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; FACILITY: MATH LIBRARY 
++ 
; ABSTRACT: 


; MTH$SSIN and MTHSCOS are functions which return the hag ed point 
; sine or cosine value of their single precision floating point argu- 
; ment (radians). The call is standard call-by-reference. 

3; MTHSSIN_R4 and MTHS$COS_R4 are special routines which are the same 
z: as MTHSSIN and MTHSCOS” except a faster non-standard JSB call is 


used with the argument in RO and no registers are saved. 
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; MTH$SINCOS is a routine which returns the floating point sine and 

; cosine value of its single precision floating point radian argument. 
; The call is standard call-by-reference. MTHSSINCOS_R5 is a special 
: routine which is the same as MTHSSINCOS, except a faster non- 
3; standard JSB call is used with the argument in RO and no registers 


; are saved. 
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; MTHSSIND and MTHSCOSD are functions which return the eating point 
; sine or cosine value of their single precision floating point argu- 


; ment (degrees). The call is standard Sptiooy rater snes. 
; MTHSSIND_R4 and MTHSCOSD_R4 are special routines which are the same 
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; as MTHSSIND and MTHSCOSD” except a faster non-standard JSB call is 
; used with the argument in RO and no registers are saved. 


; MTHSSINCOSD is| a routine which returns the floating point sine and 
3; cosine value of its single precision floating point degree argument. 
; The call is standard call-by-reference. MTHSSINCOSD_R5 is a special 
; routine which is the same as MTHSSINCOSD, except a faster non- 
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git ; standard | 488 call is used with the argument in RO and no registers 
; are saved 
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o7 : 3; VERSION: 1 
65 : HISTORY: 
6 AUTHOR : 
$00 68 ; 
i 9 ; MODIFIED BY: 
0 71: 1-1 Tryggve Fossum, 28-May-78 
000 72; o : 
0000 73 3 
44 74 : VERSION: 2 
0000 > 3 
4 £6 5 eereets 
000 77 ; AUTHOR: 
0000 78 ; BOB HANEK, 25-MAY-78: Version 2 
0000 79 ; 
0000 H+ ; MODIFIED BY: 
0000 12 
0000 82 ; Jeffrey C. Wiener 1-APR-83 
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HISTORY ; Detailed Current Edit History ~SEP-1984 MTHRTL.SRC Suvx 
4 -SBTTL HISTORY ; Detailed Current Edit History 
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ALGORITHMIC gh pes FROM tn —— 
- Uses POLYF, 8 nore accura 
2. Checks or 435 tt-] 4 before yolynomial instead of checking for 
prea er than 1.0 afterwards. 
3. Signals couplete loss of significance. 


Edit History for Version 1 of MTHSSINCOS 


Add JSB entry porate MTHSCOS_R4 and MTHSSIN_R4. Tryggve, {Uyate-76 

Change routines so as net to Use R5 ang 6. Tryggve, Jun-15-78 

Use EMODD if it > rryo9ve. ro -7 

Do CVIDF after EMODD. TNH June- 

Check argument range after EMODF. Toyaqve iy J yge-Erere 

Do CVTFD perere iy DD at CBIG. TNH 2 

Fix ov flow SIN(2**25*PI), too Large ansuers S SIN(2**1 PI) 
1=21:24. H 26-July-7 

Same fix .. ‘aeune ton arg as in 1-8. TNH ete sulyers 

Pickup arg with MOVF so reserved op check. TNH 16-Aug-78 

1 = Put version number in standard format: three digits in 
edit er intics JBS 16-NOV- 

¢ - Change, Ath s ama iates to MTHSK_SINCOSSIG. JBS 07-DEC-78 

Add o the T directive. JBS 22-DEC-78 

4 = Add RSB after Assy to MTHSSS1GNAL $9. that TAN doesn't 
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© into infinite loop. SBL 02-Feb- 
- Beclare externals. SBL 17-May-19 
- Change MTHSK_SINCOSSIG to MTHSK_SIGLOSMAT. JBS 19-SEP-1979 
- change arguaent Limit to 2**31 So as to be compatible with 
pigner Ciniere on cous tees. SBL 31-0ct-197 
tt imit go 2**30 because this routine can't handle 2**31. 
Nov-1 


1-019 - tr RINBTLAR. PACKET somehow left out! SBL 2-Nov-1979 
1-019A = Changed BLSS to BLSSU after unbiasing the exponent to exercise 
sael argument logic. This logic was never reached prior to this 
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hod fieg logic for processing reduced ar uments close to pi (sin) 
and et/ (cos) to perform the operation RARG in double 
precision. Performing the operation in single precision results 
on hag ing as much as the last 6 bits. This modification was made 
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ine ut argument with L, nitude less than 2** 
{irs word of gh $ OW OF CO entry of the €STB2 table 
from octa 5 2174 to octal 32174. This eliminated a negative 1 (sb 
iasing of the values returned by COS for small arguments. 
- RNH 12-FEB-1981 
1-20 = Added degree entry points. RNH 1-MAR-1981 
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Algor tshete differences form Version 1 
) Introduction of pm : entr onesnt 
) Elimination of t rene ction on the argument 
) Elinination of spe ATH K_SIGLOSMAT erro 
4) Introduction of a possible woe em EC. pe for radian arguments 
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; 2-004 


; ; 2-005 


pa19B6 11:58:16 EMTMRTL. SACIOVXSINCOS.MaR;1 2% 


aodt ted geass oh ARGe E §8.¢s eliminate [se bug (similar to the 
bu $ nents qoneg t 6.) RNH 1 

rg be ree SS 1Ncos. a. microvax by retsce {ne all D o. [soot ing 

pet nt tnstry ctions with equivalent G_floating point” instructions. 


Modi tied ntti & that it retucne the same values as MTH$SINcos 
in the waeerye 3 PI/4). The fix was to modify the special code 
used in — for ix! Pes that interval. The meg) i*eat ion 
of this ne. y* because there was no need to JSB t 
NEEDS bol ste -SINCOS to gon ute COS. This prob m wee \ocounred 
by nkenbeard in the ores group. JCW 
To Had the £58 error to less the cOstecl “table was changed 
from a Bais al sf degree (Hart 3850 0) to a polynosat of degree 4 
(Hart 3 The CO coefficient is weet? s CO; the earlier version used 
CO-1 to improve accuracy. Now, it is more accurate to use C0. The 
speraer sete code modifications were also made. Various comments were 
so corrected. JCW 21-Jun-84 
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Petty BECLARAT IONS = Declarative Part of Modu ets 96:0: 18 EMTHRTL. SREJOVXSINCOS MAR: 1 oe (3) 
163 -SBTTL DECLARATIONS - Declarative Part of Module | 
165 ; 
196 ; INCLUDE FILES: MTHSJACKET.MAR 
| 
168 ; EXTERNAL SYMBOLS: | 
09 190 * .DSABL GBL | 
171 -EXTRN MTHSAL_4 OV_PI 
00 178 “EXTRN MTHSSSIGNAL | 
17 ~EXTRN MTHSK_FLOUNDMAT 
000 174 “EXTRN MTHSSJACKET_TST | 
000 176 ; EQUATED SYMBOLS: | 
Henge 000 178 LF_PI_OV_4 = *X0FDB4049 : 0.78539819 | 
1D641E2 0000 179 LF-9 Pl Ov_.4 = *x 106416 : 7 06838349 
000 180; LF72"OV-PI~ = “KF 983402 : 636619772 | 
5F304004 0000 13) G.2_6v_Pi = *XCBB26DC95F 304004 = ,6366197723675813E+00 
CBE44116 0000 1 é LF_3_PT_OV_4 = XCBEG41 16 3 2.35619450 | 
53014178 0000 18 LF-5"PI“OV-4 = *X53D1417B : 3.92699075 
EDDF41AF 9000 18% LF-7~PI-OV~4 = “XEDDF4IAF : 5.49778700 | 
00004 0000 136 w_45 = *xX4334 ; First word of 45 in F-format 
00004334 0000 £187 LF_45 = *x00004334 | 
B00 0000 86188 LF_M45 = *x0000C334 
gee D65 0000 189 LF-SMALLD = *X2EE13D65 ; 180/P1*2*-12 
Ss03bBe 9000 19% Tha es = 4X88603DB6 
Nites: 0000 19¢ LF- CONVERT = anaees Spee 3 pi/180 - 24-6 | 
ZEET42ES 9000 19 [F790_0V.PI = “X2EE142E5 : 58.64789009 | 
£E10365 B208 13 LF-SMALLEST_DEG = *X2EE10365 ; 2*-128*180/PI 
om i come 
0000 198 SSFDEF : Define SF$S (stack frame) symbols 
0000 199 SPSLDEF ; Define PSL (Processor Status Longword) 
50 $ : symbols 
$444 8 3 PSECT DECLARATIONS: | 
0000000 04 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
000 05 : program section for math routines | 
0000 0? : OWN STORAGE: none | 
900 09 : CONSTANTS: | 
00 11 G_PI_Ov_2: | 
20185444 21FB4019 + 1 pags -QUAD %*X2D18544421FB4019 : 1.5707963267948966 | 
20185444 21FB4029 09 13 .s 6 aT *X2D18544421FB4029 : 3.1415926535897932 | 
21D27F 33 D97C4032 3 = # oes .GUAD © *X21D27F33097¢4032 : 4.7123889803846899 
20185444 21FB4039 a} 18 ~~ © QUAD *X2018544421FB4039 ; 6.2831853071795865 | 
| 
| 
| 
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BECLARAT IONS - Declarative Part of Modu aaa 95:98:18 ERTHRTL. SREJUVXSINCOS MAR: 1 . (5) | 
1 
§ -SBTTL COEFFICIENT TABLES - Series Coefficients 
4 | 
f 
° | 
6 3 ; Polynomial Coefficient tables for arguments in radians 
P | 
0 
6 1 COSTBR1: i COS coefficients . argunen s less than 1/2 
ce 020 § . LONG x6 FOREGS 3 C3 = -.13765 GRE-Be 
ASS5E 024 -LONG “XABA5SSE2A 3; C2 = .41664758E-01 
FFF DBFFF 028 4 ~ LONG “XFFFOBFFF : C1 = -.49999991E+00 
aes 9 0 ‘ 5 -LONG %x00004080 : CO = .10000000E+01 | 
0000004 3 § COSLENR1 = .-COSTBR1/4 | 
Ooh 8 COSTBR2: i £95 Geet ttetencs for orguaents greeter than 1/2 
836538CC 0030 39 . LONG X836538CC 3 C4 = .24379880E-04 
03C 6 0034 40 -LONG “X03C3BBB6 3; C3 = -.13886619E-02 
AASD3E2A Baee 41 - LONG Rast A 3; C2 = .41666616E-01 
€263284 003C g -LONG *%X3C263284 3: C1 = .38485437E-08 
00000000 0040 4 -LONG *%*x00000000 : CO = .00000000E+00 
00000005 ope? rt: COSLENR2 = .-COSTBR2/4 
0044 246 SINTBR: j SIN coefficients | 
BAELBAGC 0044 247 -LONG — EXBAEGBAG ; C3 = =.195066968E-03 
83023D08 0048 gag LONG %x83023D08 ; C2 = .833201595E-02 
AAAOBF2A 004C 49 -LONG “XAAAOBF2A : C1 = =. 166666508E+00 
BADEB25E 0050 250 -LONG “X8ADEB25E ; CO = -.323841887E-08 
00000004 0054 251 SINLENR = .-SINTBR/4 
0054 26 
0054 5 
0054 54 
0054 55 
0054 256 
0054 257; 
0054 258 ; Polynomial coefficients for arguments in cycles 
Bt 5 
Sad 61 COSTBC1: jcOS coefficients for arqunents less than 2/pi | 
CEAAS7OC 0054 362 LONG *XCFAAS76C : 4 = 0.35287617E-05 | 
E275 58 6 LONG onehectnas : C3 = 0°38 a7 830 bs 
EOED3D81 2f 64 -LONG “*XEOED3D81 3 C2 = 0.15854 ree 50) 
da v¢ 6 65 - LONG oa E6BF9D : C1 = -.30842513E+00 
0004080 0064 6 -LONG %*x00004080 : CO = 0.10000000E+01 
0000000 be tf COSLENC1 = .-COSTBC1/4 
06 $9 COSTBC2: iCOS coefficients for ar yn nts greater than 2/pi 
CFAA376C 006 0 .LONG  *XCFAASTEC ; ch o 8. 87617 “65 | 
e673 of 71 . LONG one RZ RAMS 3; C3 =e. ¢ 79306 -85 
E +) p81 7 -LONG “XE 081 3; C2 = 0.15854 ¢ E-01 
4F 6F 74 7 -LONG “X4F 52BE6F 3; Cl = 9° Becher 6E-01 
008 3 78 74 LONG mS 0000 : CO = 0.00000000E+00 
00000 As 75 COSLENC2 = .-COSTBC2/4 
o7¢ 77 SINTBC: : SIN coef for arg in cycles 
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FUNCTIONAL DESCRIPTION: 


The SIN, COS and SINCOS routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = 11e(2*pi) + Ie(pi/4) + YI, 


0 =< I < 8 and 0 =< Y1 < pi/4. Since SIN and 
follows that 
/4) 
/ 


t 
pi/4) + ¥1)) and 
pi/4) + Y1)). 


where I1 and I are integer 


ers 
COS have a period of 2*pi i 
SIN(x) = SIN(I#( 
COS(x) = COS(1*¢ 


Using the trigonometric identities for the sum and difference of two angles, 
the following table can be generated: 


If I = then SIN(x) = and COS(x) = 


-SINCY1) 
-COS(pi/4-¥1) 
-COs(Y1) 
4-Y1) “Sisrerert 
SINCY1) 
4-Y1) COS (pi/4-Y1) 
Let Y be defined as Y = Y1 is even and Y = pi/4 - Y1, if I is odd, then 
each entry of the above table is of the for +/-SIN(Y) or +/-COS(Y). Based 
on the above remarks, the SIN, COS and SINCOS routines process the input 


argument x, to obtain I and Y, and based on | selects a suitable polynomial 
approximation, p(Y), to evaluate the desired fuction. 
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INPUT PARAMETERS: 
LONG 
x 


sine 
cosine 


45419494 
0 
0 


1*LONG 3; x is input angle in radians 
: sine is SIN(x 
*LONG 3 cosine is COS(x) 


nana 
wns 
7 
rm 
o 
z 
a 


SOOOoCOoOoOoOOooOooOoOoooooooooooooOoOoOOOOOOOOOODOOOOOOoOO 
SOGCOGOCOCOSOOSOSSOOCOSCSOSOSOSOSCSOOOSOSOSOSOOOOOO OOOO OOOCOOOOOSOOO 
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| 
| 
| 
| 
Cees ioe : ~ COs(y) 
i/4-¥1) SIN(pi/4=Y1) 
| 
| 
| 
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XSSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 21 AX/VMS Macro v04-00 
yyaO8 ATHSSINCOS = = Radian arguments 6-SEP=1 1984 tf 98: 18 UMTHRTL. SREJUVXSINCOS MAR: 1 
C 61 
C 86 
C 65 ; 
: rt: ; Return sine and cosine of argument 
C 66 * 
gh 
003C 5 oe ENTRY MTHSSINCOS, “M<R2, R3, RS, RSD 
Be % MTHSFLAG_JACKET 
6D 00000000 ' GF 9E OC MOVAB G*MTHSSJACKET_HND, (FP) 
0c ; set handler address to jacket 
0c9 ; handler 
00C9 
Bate 71 
50 04 6C 50 C9 MOVF ax(AP), RO 
oppoonse"et 16 00C 7 JSB MTHSSINCOS_RS 
08 BC 0 pO 00D 74 MOVL RO, a@sine(AP) 
0c BC 51 £4200 OOD 75 MOVL R1, a@cosine(AP) 
04 00DB 76 RET 
00DC 77 
00DC 78 
00DC 79 
00DC 80 -SBTTL MTHSSIN 
00DC 81 
00DC +f 3 
00DC 83 ; Return sine of argument 
00DC 84 ; 
00DC 85 
00DC 86 
001¢ OODC 387 eENTRY MTHSSIN, “M<R2, R3, R4> 
OODE $65 
BODE 89 MTHSFLAG_JACKET 
6D 00000000 ' GF 9E OODE MOVAB G*MTHSSJACKET_HND, (FP) ’ 
00E5 3 set handler address to jacket 
O0E5 3; handler 
O00E5 
OO0E5 90 
50 04 B8C 50 OO0E5S 91 MOVF ax(AP), RO 
00000221°EF 16 Ope e 9 JSB MTHSSIN_R4 
04 OEF 9 RET 
oro 94 
f 95 
OF 38 
OF 9 -SBTTL MTHSCOS 
OF 98 
+4 99 ; ; 
OF 400 ; Return cosine of argument 
Bare 401 ; 
OF £86 
of 40 
001C : i0s -ENTRY MTHSCOS, “M<R2, R3, R4> 
a 406 MTHSFLAG_JACKET 


G6 2 
XSSINCOS ; Floating Point Sine, Cosi nd Si 16-$ 4 21 AX/VM v04-00 P 1 
peat Tus eos "s Ps ee ae nee g-8EF =138% 9¢: 98: 18 cya yie SUVX XSINCOS.MAR; 1 - (8) 
6D 00000000 ' GF 9E F MOVAB G*MTHSSJACKET_HND, (FP) 
F 3; set handler address to jacket 
: ; handler 
F 407 
50 04 BC 50 F9 $38 MOVF ax (AP) 
000002D6"EF 16 oF 409 JSB MTHSCOS_ re 
04 18 410 RET 
01046 411 
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00000008 
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MTHRTL.SRCJUVXSINCOS .MAR; 1 
-SBTTL MTHSSINCOSD Degrees 


FUNCTIONAL DESCRIPTION: 


The SIND, COSD and SINCOSD routines are based on octant reduction. Given an 
argument, x, it is written in the form 


x = 11*#360 + 1*45 + ¥1, 


where I1 and I are integers =< 1 < 8 and 0 =< Y1 < 45. Since SIND and 
cOSD have a period of it Sits that 


SIND(x) = SIND(1*45 + Y1) and 
COSD(x) = COSD(1*45 + 1). 


oonny the trigonometric identities for the sum and difference of two angles, 
) 


the following table can be generated: 
If l= then SIND (x) = and COSD(x) = 
0 SIND(Y1) COSD(Y 
1 COSD(45-Y1) SIND (43=¥1) 
¢ Cosp(Y1) ~SINDCYI 
SIND (45-¥1) ~tosptasev1) 
4 -SIND(Y1) -COSD(Y1) 
5 -C0SD(45-Y1) -SIND(45-Y1) 
6 -COSD(Y1) SIND(Y1) 
7 -SIND(45-Y1) COS (45-1) 
Let Y be defined as Y = ¥1 if I is even and Y = 45 - if I is odd, then 


each entry of the above table is of the for S/=SINCY)” pf +/- -Costy). Based 

on the above remarks, the SIND, COSD and SINCOSD routines process the in put 

argument ee to obtain I and Y, and based on I selects a suitable eel ynental 
oO 


approximation, p(Y), to evaluate the desired fuction. 
LONG = 4 
sind = Wei ee 
cosd = 5*LONG 


2 


t+ lt 


Q0000000'GF 9E 


50 Oe BC 50 
QOOOOS5E*EF 1 
He 
C BC D 
04 
001C 


Q0000000'GF 9 


50 04 BC 
QOOOOSAE "EF 


00000000'GF 9E 
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NCO D = Degrees 


Point Sine, Cosine ond thous i ob th 9 OG: 96: 1 ae Be 


6-SEP=19 MTHRTL 


sENTRY MTHSSINCOSD “M<R2, R3, R4, RS> 


; set handler address to jacket 


MTHSFLAG_JACKET 
MOVAB G*MTHSSJACKET_HND, (FP) 
; handler 
MOVF axa RO 
JSB aTNssiNcosD R5 
MOVL RO, a@sind(AP) 
MOVL R1, @cosd(AP) 


.ENTRY MTHSSIND 
MTHSFLAG_JACKET 
MOVAB G*MTHSSJACKET_HND, (FP) 


3; handler 


“M<R2, R3, R4> 


; set handler address to jacket 


MOVF § aXx(AP), RO 
JSB MTHSSIND_R4 


RET 


.ENTRY MTHSCOSD 
MTHSFLAG_JACKET 


“M<R2, R3, R4> 


; set handler address to jacket 


MOVAB G*MTHSSJACKET_HND, (FP) 
3; handler 
MOVF aX(A 
J$8 minscogpe Fi 
RET 


v04-00 Page \é 
juvx SINCOS.MAR; 1 (10) 
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J.2 
e608 ® frnsSTNcos_ A "t S1P** Cosine and Sincos 1p s En 1986 17:09:16 EMTMRTL. SRCIUVASINCOS.maR;1 =o 


1? $30 -SBTTL MTHSSINCOS_RS 
14 138 ; This routine computes the SIN and COS of the F-format value of RO. The 
14 495 ; computation is performed one of three ways depending on the size of the 
1? t3e ; input argument, X: 
12 138 ; 1) it i < pi/4, then X is used directly in polynomial approximation 
: ) an ° 
14 $38 5 2) If pi/4 =< ixt < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
14 499 ; to reduce the argument to an equivalent argument in radians, Y, and 
14 0; the octant, I eantetanne the argument. Y is then evaluated in two 
14 1g go" znoe als chosen as a function of I, to compute SIN(X) and COS(X). 
14 § : 3) If Sepi/4 =< iX!, then the subroutine REDUCE_LARGE is called to 
14 3 reduce the argument to an equivalent argument in cycles, Y, and the 
14 04 ; octant, I, contanining the argument. Y is then evaluated in two 
12 82 5 polynomials chosen as a function of I, to compute SIN(X) and COS(X). 
14 09 MTHSSINCOS_R5:: 
54 50 50 Bie 08 MOUF RG 7 R4 = X 
OF 18 0148 509 BGEQ POS SINCOS ; 
00000157°EF 1 8188 219 JSB INTOS ; RO = SINCIX!), R1 = COS(X) 
~~ ww § 15 11 MNEGF RO, RO ; RO = SIN(X) 
0 Oise 818 zy 
0157 514 SINCOS: 
50 8000 8F AA 0157 515 1CW #*x8000, RO ; RO = IX! 
015C¢ 218 POS_SINCOS: 
50 OFDB4049 BF 51 B138 51 CMPF #LF_PI_OV_4, RO : Compare pi/4 with :X! 
2F 614 «(016 518 BGTR SMACL_SINTOS 3; No argument reduction is necessary 
50 31D641E2 8F 51 0165 233 CMPF #LF_9_PI_OV_4, RO ; Compare 9*pi/4 with ix! 
03 18 016C st BGEQ = 
oosc =O «5311 Bigs 521 BRw LARGE_SINCOS 3; Use special logic for iX! > 9*pi/4 
171 : 5 ; 
4 : : 3 pi/4 =< IX! < 9epis4 
O000044D'EF 16 8171 258 is: JSB REDUCE _MEDIUM ; Medium argument reduction routine 
0177 52 3; R3/R4 = Y = reduced argument 
0177 8 ; R2 = octant 
7E 53 7D 0177 9 MOVa R3, -(SP) 3; Save reduced argument on stack 
52 ODD OI7A 0 PUSHL R2 ; Save octant bits on stack 
OOOOO2FS*EF 16 O17C 1 JSB M_COS 3; RO = COS(X) 
é E 0 HT 3 MOVL (SP)+, a : Re = Octant bits 
3 1 mova (SP)+, R ; R5/R4 = reduced argument 
DD 0188 4 PUSHL RO 3; Save COS(X) on stack 
OOOO024C'EF 16 O18A 5 JSB M_SIN ; RO = SIN(X) 
51 E B? 139 § aoe (SP)+, R1 3; R1 = COS(X) 
194 3 ; : 
3; Logic for small arguments. iX! < pi/4. 
13 dp Logic f Ll t iMt < pi/4 
194 541 ° 
194 4g SMALL_SINCOS: 
50 4000 8F B1 0194 4 CMPW #*x4000, RO : Compare 1/2 with ix! 
A 19 Q199 4G BLSS $ ; Sufficent overghang not available 
50 3A80 BF 861 18 45 CMPW = #*X3A80, RO ; Compare with 2°-12 
1F «618 «(OA 46 BGEQ 1$ ; No polynomial evaluation is needed 
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acro Vv04-00 
- MTHRTL.SRCJUVXSINCOS.MAR; 1 
OVL. RO, R4 R4 = IX! 
MULF3 RG. R4, RS + RS = Kex 
POLYF RS. #COSLENR1-1, COSTBR1: RC = COS(X) 
PUSHL 3; Save COS(X) on stack 
POLYF R5, #SINLENR-1, SINTBR : RO = g(xe2) 
MULF R4. R + RO = #g(X°2) 
ADDF R4.R : RO = SIN(X) 
MOVL (SP)+, R1 + R1 = COS(X) 
RSB 
RSB 
PUSHL RO » Save ix! on stack 
Replace 
CVTFD RO, R3 
MULD2 R3, R3 
CVTFG RO, R4 : RG/RS = IX! 
MULG2 R4, RG + RG/R5 = K*2 


R3/R4 is needed in its G_floating Copregentet ten for NEEDS_DOUBLE_SINCOS. 
The F Fpoet hae /eorenent ation of R3/R4 is also needed after the return 
from NEEDS _DOUBLE_SINCOS. The following code will accomplish these goals. 
ihe only ting of code that is being replaced in the D_ to G_ conversion 
sa . 


Note, the BICL3 used below is necessary since a CVIGF yields a rounded 
value instead of a desired truncated value is obtained. The rounding bit 
(bit 12 of the second longword) plus all lower order bits of word 3, that 
is, the first word of the second longword are zeroed. 

MOVL R4, R3 


BICL3 #*XOOOO1FFF, R5, R4 


R3 = first longword of G_floating 

representation of x* 

Rounding bit and all lower order 
bits are zeroed. Original second 
longword of G_floasting rep of 
X*2 is saved Tn RS 

Save X*2 on the stack 

R3/R4 is G_floating rep of x*2 


CVIGF R3, =(SP) 
R5, R4 


MOVL 
PUSHL R353 3; Save x*2 on stack 
JSB artes DOUBLE _SINCOS ; Use special logic to obtain overhang 
MOVL RO, R5 3; Save COS(X) in RS 
MOVL (SP)+, RO > RO = x42 
POLYF RO, #SINLENR-1, SINTBR : RO = : X*2) 
MULF (SP), RO : RO = #g(X°2) 
ADDF (SP)+, RO 3 RO = SIN(X) 
MOVL 5, Rf > R1 = COS(X) 
RSB 
LARGE_SINCOS: . 
JSB REDUCE _LARGE ; R3/R4 = reduced argument (in cycles) 


. 2 
YESS INCOS ; Floa sting ie Sine, Cosine and Sincos 16-SEP-1984 03:06:12 AX/VMS Macro V04-00 P 
| 2-00 MTHSSINCO 6-SEP=1984 11:29:16 CMTHRTL.SRCIJUVXSINCOS.MAR; 1 
0 oe 3; R2 = octant bits 
2¢ 4 0 605 PUSHL § 3; Save octant bits on stack 
7 D 0 on6 MOVa R353, -(SP) 3 Save reduced argument on stack 
00000 32'6F 16 08 607 JSB L_¢os ; RO = COS(X) 
—E 7D E one MOVQ (SP)+, R3 ; . duced argument in R3/R4 
dO 11 + 60 MOVL (SP), R2 3 = octant bits 
gf dO 14 610 MOVL RO, (SP) ; RS = COS(X) 
00000 AA'Er 16 0217 611 JSB L_§INn : RO = SIN(X) 
51 E ? 1D ol¢ MOVL (SP)+, R1 ; R1 = COS(X) 
20=—s «61 RSB 
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$ Page 1 
6-SEP-1 MTHRTL.SRCJUVXSINCOS .MAR; 1 ° H$) 
-SBTTL MTHSSIN_R4 


a id 


ro 


This routine computes the SIN of the F-format value of RO. The computation 
Hs performed one of three ways depending on the size of the input argument, 


1) If ix! < pi/4&, then X is used directly in a polynomial approximation 


fr) ° 

2) If pi/4 =< ix!t < 9epi/4, then the subroutine REDUCE_MEDIUM is called 
to reduce the argument to an equivalent argument in radians, Y, and 

the octant, I, containing the argument. Y is then evaluated in a 

pg 'xnony et chosen as a function of I to compute SIN(X). 

f 9epi/4 =< iX!, then the subroutine REDUCE_LARGE is called to 
reduce the argument to an equivalent argument in cycles, Y, and the 
octant, I, contanining the argument. Y is then evaluated in a 
polynomial chosen as a function of I to compute SIN(X). 


MTHSSIN_R4:: 
TS 


3) 


50 53 TF RO ; Check the sign of RO 
OF 18 BGEQ POS_SIN ; 
0000022F ‘EF 16 JSB IN 3; RO = SINCGIX!) 
50 50 3 MNEGF RO, RO ¢ RO = SIN(X) 
0 RSB 
SIN: 
50 8000 8F BICW #*x8000, RO ; RO = ix 

POS_SIN: 


"CMPF «= ss #LF_PI_OV_4, RO ; Compare pi/4 with !Xx! 
BGTR SMACL_S1 ; No argument reduction is necessary 
CMPF #LF_9_PI OV_4, RO ; Compare 9*pi/4 with {Xi 


AA 

50 OFDB4049 8F i 
1 

y BLSS LARGE-SIN ; Use special logic for ‘Xi > 9*pi/4 


$e 
50 31D641E2 8F 
5E 


: pi/4 =< IX! < 9epis4 
JSB REDUCE_MEDIUM 


0000044D°EF 16 


PDA-_DPLAAAAAAAAAAA AAAS AAA AA AAA AAA AAA AAA AAAAAAAAO ASO AO 
DPV SVT EB IS BS ES SE EE NANI WII NIAAA ROPPIPNIPUNYINUN 
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FUFUSUSISISIN E™ BE EE ES EE ENA DINPIPUPIPOPIPIPIPIPCPIPIPIPININPIPIPIFININYD st 
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; Medium argument reduction routine 
s ay 7 ; reduced argument 
3 = octan 
07 OO 52 8 M_SIN: CASEB R2, #0, #7 ; Branch to one of four polynomial 
3 evaluations depending on the 
0586' 1$: -WORD P_SIN_R-1$ 2 octant bits. 
0514° «WORD P_COS_R-1$ 
0514' «WORD P_COS_R-1$ 
581° .WORD NZSINZR-1$ 
5B1' -WORD N7SINR-1$ 
360° .WORD N7COS"R-1$ 
560° 2e 661 «WORD N_COS_R-1$ 
586° 66 - WO P_SIN_R-1$ 
6 66 
6 664 ; : , 
6 665 ; Logic for small arguments. iX! < pi/4. 
6 066 : 
6 66 
6 668 SMALL_SIN: 
50 4000 8F 81 6 06 CMPW = #*X4000, RO ; Compare 1/2 with {X! ‘ 
A 19 6 670 BLSS 2s ; Sufficent overghang not available 
50 3A80 BF BI 6 671 CMPW #*x3A80, RO 3; Compare with 2-12 


N 2 
VXSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:06:1 AX/VMS Macro V04-00 Pa 17 Vv 
pet ATHSSIN_RE pees 7 96:96:18 EMTHRTL. SRETOVXSINCOS MAR: 1 - (12) Ss 
1 is 6C ore BGEQ 1$ 3; No polynomial evaluation is needed 
26 1) of of MOVL Re. RG 3 RO = Xx 
44 74 MULF RO,RO ; RO = Mex | 
FDCA CF 2 ? 2? oe of? fore aye “pets SINTBR ; RO = q(x*2) 
5 4 40 7D 699 ADDF R4, R ' RO = SIN(X) 
0 ? ors 1$: RSB 
50 DD 1 o80 2$: PUSHL RO 3 Save ix! on stack 
28 50 99FD ; 681 CVTF RO, R : RO/R1 = IX! 
0 50 44FD ? ¢ MULG RO, R ; RO/R1 = x*2 
B 684 ; Note, the BICL3 used below is necessary since a CVTGF yields a rounded 
B 685 ; value instead of a desired truncated value is obtained. The rounding bit 
B 686 ; (bit 12 of the second longword) plus all lower order bits of word 3, that 
oe ets ; is, the first word of the second longword are zeroed. 
8B 689 ; RO = first Longword of G_floating FF 
OSe0 690 : representation of X*2 
51 51 QOOO1FFF 8F CB 8 691 BICLS #*XOOOOIFFF, R1, R1 ; Rounding bit and all lower order 
029 $26 ; bits are zeroed. Original second 
029 69 ; longword of G_floasting rep of 
029 694 ; &*2 is saved in R5 
50 50 33FD 029 695 CVTGF RO, RO 3; RO = x* 
FDA7 CF 03 50 55 0297 696 POLYF RO, #SINLENR-1, SINTBR ; RO = Q(X*2) 
50 6€ 4&4 029d 697 MULF2 (SP), RO + RO = X*qQ(X*2) 
50 BE 40 OSAg 698 ADDF (SP)+, RO 3: RO = SIN(X) 
0 Ose 900 chee 
BSA 701 LARGE_SIN: 
OOOOO4C1'EF 16 rd re JSB REDUCE_LARGE ; een 8 Peeters erpmem (in cycles) 
3 = oc 
53 DS O2AA 704 L_SIN: TSTL R3 ; Check for degenerate case 
14 13 3 ry ree BEQL DEGENERATE_CASE_SIN 
07 OO 52 &8F QO Y a4 CASEB R2, #0, #7 
43h B 709 1$: «WORD P_SIN_C-1$ 
0578" 0284 710 -WORD P7COS"C-1$ 
44 os 711 -WORD P_COS_C-1$ 
619° 028 re -WORD P_SIN_C-1$ 
614° O2BA 71 «WORD N_SIN_C-1$ 
05C5' O2BC 714 -WORD N7COS"C-1$ 
B20? 8 BF 715 «WORD N_COS_C-1$ 
614" 02C ig WORD N7SIN7C-1$ 
O55 718 
: At, DEGENERATE_CASE_SIN: 
2 ga C 721 BICB #i, R2 3; Compute index as (R2 - 1)/2 
52 FF af c 02cS 7 : ROTL #1, R2, R2 | 
03 +00 ae cA f ; ASEB R2, #0, #3 
718" OSCE 795 1$:  .WORD-P_ONE-1$ | 
724° D 7 § -WORD URFL -1$ 
71C° D 7 -WORD WN _ONE-1$ 
724° D4 728 -WORD URFL -1$ 
| 
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IN_R 66-SEP-1984 MTHRTL.SRCJUVXSINCOS MAR; 1 
1 
$e 
D6 f : .SBTTL MTHSCOS_R4 
D6 «735 ; This routine computes the COS of the F-format value of RO. The computation 
D660 7 § ; is performed one of three ways depending on the size of the input argument, 
4 4 P ; X. The processing is the same as described for MTHS$SIN_R4. 
oe 739 
5053 Ospe Fey NSCOS MRE ko Check f d d 
: Check for reseved operan 
50.8000 BF AA D8 748 Bicw #*X8000, RO SRO IX: F 
50 OFDB4049 BF 51 ODD 74 CMPF = #LF_PI_Ov_4, RO > Compare pi/4 with {Xx! 
$3 14 O2E4 744 BGTR SMACL_COS ; No argument reduction is necessary 
50 31D641E2 BF 51 O2E6 745 CMPF  #LF_9"PI_OV_4, RO : Compare 9*pi/4 with Xx! 
30 019 8 ED eg BLSS LARGE-COS 3 Use special logic for ix! > 9epi/4 
O2EF 748 ; 
8 4 rep 3 piss =< IK! < 9episs 
0000044D'EF 16 O2EF 751 JSB REDUCE _MEDIUM ; Medium argument reduction routine 
O2F5 r3¢ ; R3/R4 = Y = reduced argument 
O2F5 75 ; R2 = octant 
07 OO 52 8F O2FS 754 M_COS: CASEB R2, #0, #7 ; Branch to one of four polynomial 
O2F9 755 $ evaluations depending on the 
46B' O2F9 756 1$ -WORD P_COS_R-1$ : octant bits. 
0508° OeFB 757 .WORD N7SIN-R-1$ 
0508' O2FD 758 «WORD N_SIN_R-1$ 
0487' ogee 759 «WORD N_COS_R-1$ 
0487" 0301 760 -WORD NZ COS-R-1$ 
050D" 0303 761 «WORD P_SIN_R-1$ 
050D' 0305 76 «WORD P_SIN_R-1$ 
0468' 0307 76 .WORD PTCOS"R-1$ 
0309 764 
0309 765; ‘ 
0309 766 ; Logic for small arguments. iXi < pi/4. 
0309 767 ; 
0309 768 
Oab8 769 SMALL_COS: : 
50 4000 8F B61 0309 770 CMP #*x4000, RO ; Compare 1/2 with iXi_ ; 
07 ~=«+18 O30 44 BGEQ 1% : Sufficent overghang is available 
031 ae: : CVTFD RO, R3 > R3/R4 = IX! 
$318 774 ; 
53 50 99FD 1 775 CVTFG RO, R3 3 R3/R4 = IX! - 
045D 31 0314 776 BR NEEDS DOUBLE ; Use special logic to obtain overhang 
50 3a80 BF 861 17 777 1$: CMPW #*X3AB80, RO ; Compare with 2*-1 : 
A 18 1C¢ ore BGEQ : No polynomial evaluation is needed 
0 0 44 1. BE ks MULF RO,RO 3; RO = Xx 
FCF9 CF 0 oe 1 780 POL YF RO, #COSLENR1-1, COSTBR1; RO = COS(X) 
05 7 , 1 RSB 
50 08 50 ; 7 ; 2$: MOVF #1.0, RO ; RO = 1.0 = COS(X) 
05 B 784 RSB 
C 785 
ie : $ 
C 787 LARGE_COS: 
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QOOOO4C1"EF 16 If $ JSB REDUCE _LARGE ; R3/R4 = reduced argument (in cycles) 
; R2 = octant bits 
53 rf ? LCOS: TSTL R3 Check for degenerate case 
14~—COYT 1 a BEQL DEGENERATE_CASE_COS 
07 OO 52 &F 6 £38 CASEB R2, #0, #7 
4F3' O33A 794 18: .wORD P_Cos_C-1$ 
591° Cc 795 WORD P_SIN_C-1$ 
58c' 196 “WORD NTSIN7C-1$ 
530° 4 79 «WORD N_COS_C-1$ 
053p" 034 198 “WORD N=COS"C-1$ 
O58¢: 44 79 “WORD NZSIN7C-1$ 
591° 6 00 «WORD P_SIN-C-1$ 
O4F3" 034 01 “WORD P7COS~C-1$ 
034A o 
034A 0 
034A 04 DEGENERATE_CASE_COS: 
034A Bh? 
52 BA 034A 06 BICB #1 ; Compute index as (R2 = 1)/2 
52 52 FF BF 9C 0340 80 ROTL #1, Ree, R2 
03 00 8F 0352 808 CASEB R2, #0, #3 
0356 809 
069C" ba28 si? 1$: -WORD UNFL -1$ 
0694" 035 11 -WORD WN_ONE-1$ 
069C* O35A aig WORD URFL -1$ 
0690" 035c 81 [WORD P_ONE-1$ 
O35E 814 
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} 16 -SBTTL MTHSSINCOSD_RS 

E is ; This routine computes the SIND and COSD of the F-format value of RO. The 

E 19 ; computation is performed one of two ways depending on the size of the input 

; Y ; argument, X: 

E ¢ 3 1) If iki < 65, then X is used directly in polynomial approximation 

E 3 of SIND and COSD. 

E 4; 2) If 45 =< ixi, then the subroutine REDUCE_DEGREES is called to reduce 
5E ee the argument to an equivalent argument in degrees, Y, and the 
5E § 3 octant, I, containing the argument. Y is then evaluated in two 
ef 3 polynomials chosen as a function of I, to compute SIND(X) and 

E 8 ; OSD(x). 

O35E 4 
O35E Q MTHSSINCOSD_RS:: 
50 6553 8 43 1 TSTF RO 3 
OF 18 6 § BGEQ SINCOSD ; 
50 8000 8F AA OQ os ; BICW #*x8000, RO 3: RO = IX! 
00009 Ter 1 036 4 JSB SINCOSD 3; RO/R1 = SIND( IX!) /COSD(iX!) 
50 0 52 036d 35 MNEGF RO, RO ; RO = -SIND(IX!) 
05 0370 o3$ RSB 
0371 37 
0371 838 SINCOSD: 
50 4334 8F 61 0371 839 CMPW #W_-5, RO ; Compare 45 to !X: 
1€ 14 0376 840 BGTR SMALL_SINCOSD : special processing for small arg 
OOOOO6E2"EF 16 0378 841 JSB mgouce DEGREES ; R3/R4 = octant/reduced argument 
7E 3 7D O37E eee mMOva R3, -(SP) 3; Save octant bits and reduced arg 
00000426'EF 16 0 84 JSB EVAL_COSD ; RO = COSD(Y) 
as. 8 fF 844 mova (SP), R3 ; R3/R4 = octant/reduced argument 
50 Dd 0 845 PUSHL RO ; Save COSD(Y) 
OOOOOSCE*EF 16 O 846 JSB EVAL_SIND ; RO = SIND(Y) 
51 8E dO 0 847 MOVL (SP)#, R1 3; R1 = COSD(Y) 
05 0 848 RSB 
0 849 
0 850 
0 851 SMALL_SINCOSD: 
50 0 S26 PUSHL RO ; Save argument 
000004 3A'EF 0 8 JSB SMALL_COSD : RO = COSD(:X!) 
55 50 0 54 MOVL 3: RS = COSDCIX!) 
0 855 MOVL  (SP)+, RO ; RO = argument 
0 56 JSB SMALL_ SIND 3; RO = SIND(X) 
0 57 MOVL R5, RT : R1 = COSDCIX!) 
0 58 RSB 
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.SBTTL MTHSSIND_R4 


16 CMTHRTL.SRCJUVXSINCOS.MAR; 1 


; This routine computes the SIND of the F-format value of RO. The computation 
; is performed one of two ways depending on the size of the input argument, X 


3 1) if ist < 45, then X is used directly in polynomial approximation 
: fr) . 
; 2) If 45 =< ixt, then the subroutine REDUCE_DEGREES is called to reduce 
; the argument to an equivalent argument in degrees, Y, and the 
: containing the argument. Y is then evaluated in two 
chosen as a function of I, to compute SIND(X). 
MTHSSIND_R4:: 
TSTF ; RO = X 
BGEQ : 
JS 3 RO = SIND(IX!) 
MN 3 RO = -SIND(IX:) 
RS 
NEG_SIND: 
BIC ; RO = IX! 
POS_SIND: 
CMPW #w_45, ; Compare 45 to {X! 
BGTR SMALL_SIND : sgec tat processing for small arg 
JSB REDUCE_DEGREES ; R3/R4 = octant/reduced argument 
EVAL_SIND: 
CASEB R3, #0, #7 
1$: .WORD P_SIN_6-1$ 
eWORD P_COS_D-1$ 
«WORD P_COS_D-1$ 
.WORD P7SIN7D-1$ 
.WOR N_SIN_D-1$ 
«WORD N_COS_D-1$ 
- WOR N_COS_D-1$ 
i N_SIN_D-1$ 
SMALL_SIND: 4 
CMPF #LF_SMALLD, RO ; Compare 180/pi*2*-12 with ix: 
BGTR 1$ 3; No polynomial evaluation is 
MOVL 3 necessary 
BRW : 
1$: TSTF ; Check for zero 
4 BEQL 3$ ; Return if RO = 
CMPF #LF_SMALLEST_DEG, RO ; Check for possible underflow on 
BLEQ 2$ ; conversion to radians ; 
BRW UNFL ; Underflow will occur on conversion 
2s MULF3 #LF_CONVERT, RO, R1 ; RI = teilize ~ 2°=6) tix! 
SUBW x300 3 RO = iXit2*- 
aoe 3: RO = SINDCIX!) = (pi/180)ix: 
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6-SEP=-19 
-SBTTL MTHSCOSD_R4 
; This routine computes the COSD of the F-format value of RO. The eer 
: is performed one of two ways depending on the size of the input argument, xX 
; Details are given in the d scussion on MTHSSIND_R4. 


MTHSCOSD Beit 


MTHRTL.SRCJUVXSINCOS.MAR; 1 


TSTF R ; Check for reserved operand 
BICW #*x8000, RO ; RO = 
CMPW 45, & ; esata ‘45 to IX! 
BGTR SMALL_COSD ; 
JSB REDUCE _DEGREES ; R3/R4 = octant/reduced argument 
EVAL_COSD: 
CASEB R3, #0, #7 
1$ .WORD P_COS_6-1$ 
“WORD PTSIN“D-1$ 
[WORD NTSIN“D-1$ 
“WORD N7COS"D-1$ 
“WORD N7COS~D-1$ 
-WORD N_SIN_D-1$ 
“WORD P7SIN7D-1$ 
«WORD P_COS_D-1$ 
SMALL_COSD: 
CMPF #LF_SMALLD, RO : Compare 180/pit2*=12 with ix! 
BGTR 1$ 3; Check if polyinomial evaluation is 
MOVL RO, RG r necessary. 
BRW P_COS_D 3 POLY neede d 
1$: MOVF aT, RO ; RO = 1, = COSD(IX!) 


] 
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Wu08 AtHscosp. RG paar 7 94:98:18 ERTHRTL . SRETUVXSINCOS.MAR: 1 . AS). 
44D v4 
44D t 
44D 94 -SBTTL REDUCE_MEDIUM 
re 
440 9 4 ; This routine assumes that the absolute value of the argument, X, is in RO 
46D 955 ; and that pi/4 =< ix! < 9epi/4. It returns a double precision reduced 
ret gee 3; argument, Y, in R3/R4, and the octant bits in R2. 
44D 32$ 3; The reduced orgenent is obtained by locating the octant that X is in through 
rt 82 3 a@ binary search and then subtracting off a suitable multiple of pi/2 | 
tcp 989 | 
44D 960 REDUCE_MEDIUM: 
0440 961 ; 
44D 366 ; Omit the Line of code: 
44D 9635 ; CLRL a 3 RO/R1 = IX! | 
44D 964 ; 
50 53014178 ef 51 0440 965 CMPF #LF_5_PI_OV_4, RO 3 
3 615 #0454 966 BLEQ 5$ 3 IM! >= Sepi/4 
50 CBE44116 BF 51 0456 967 CMPF #LF_3_PI_OV_4, RO 
wm: Beee 43 BLEQ 3$ 3 IX! >= 3epi/4 
Beer a : SUBD3 D_PI_OV_2, RO, R3 z R3B/R4 = IX! = pi/2 | 
50 50 99FD O045F 443 CVTF RO, RO 3 RO/R1 = {X! | 
53 50 FB9B CF 43FD 0463 97 SUBG3 G_PI_OV_2, RO, R3 : R3/R4 = IK! = pi/2 
06 18 O46A 974 BGEQ 2s 3 | 
52. 01 o? 046C 975 MOVL #1, R2 3; Octant bits = 001 
05 O46F 976 RSB 
470 977 | 
52 02 00 047 978 2$: MOVL #2, R2 ; Octant bits = 010 
05 047 979 RSB 
474 980 | 
474 981 ; old SUBD3 D_P1,RO,R3 
0474 oa8 : 
0474 98 
50 50 99FD 0474 984 3$: CVTF RO, RO : RO/R1 = {Xi | 
53 50 FB8B CF 43FD 0478 985 SUBG G_P1I, RO, R3 3 R3/R4 = IX! = pi 
04 18 O47F 986 BGEQ 4 3 2 
52 03 » 481 987 MOVL #3, R2 3 Octant bits = 011 
an 
52 04 0 8 H 990 4$: MOVL #4, R2 ; Octant bits = 100 
" : 44 RSB 
50 EDDFSIAF BF 51 0489 993 5$: CMPF #LF_7_PI_OV_4, RO 3 
360COYS rs ope BLEQ 7$ 3 Kt >= 7epi/4 
rh p44) ; SUBD3 D_3_PI_OV_2, RO, R3 3 R3B/R4 = IX! = 3epi/2 
50 50 99 9 998 CVTF RO, RO : RO/R1 = IX : 
53 50 FB7S CF 43 496 999 SUBG3 6G_$_PI_OV_2, RO, R3 > R3/R4 = IK! = 3epi/2 


Octant bits = 101 


FD 
FD 1 
4 1 49D 1 o BGEQ 6$ 
52 5 3 49F : § MOVL #5, R2 
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MOVL #6, R2 ; Octant bits = 110 
RSB 


atin 
E_ME 
1004 
tA 1908 
ag 1 
ra i ; old 
4A? | 03 : SUBD3 D_2_P1, RO, R3 : RB/R4 = IX! = 2epi 
50 50 99FD O4A7 1011 7S: CVTFG RO, RO ; RO/RI = IX! 
5350 £868 CF 43FD 04 1018 SUBG3 GG Pl. RO, R3 + R3/R4 = IX! = 2spi 
5 07 ~~»v0 439 + 9 MOVL #7,°R 3; Octant bits = 11 
OF tee 101s oie 
52 DS 0486 1046 8$: CLRL Re ; Octant bits = 000 
05 0488 101 RSB 
0489 1018 
4B9 101 
4B9 1920 
0489 1021 
0489 1022 
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wyaO8 REDUCE LARGE Fe obe 94: 9: 18 YRTARIL SRC JUVXSINCOS.MAR; 1 . 38), 
4B9 1024 .SBTTL REDUCE_LARGE | 
a | 
489 1 ; ; This routine is used to reduce Large arguments (iX! >= 9*pi/4) modulo pi/4. 
4B9 1 8 ; It returns the reduced ecqument. ¥ in R35/R4 in units of cycles, and returns 
ro } 4 ; the octant bits, I, in R2. 
4B9 1; The method of reduction is as follows: | 
489 1 § ; x#(4/pi) = s Sr yt eed 3H where n an integer and 1/2 =< f <1 
489 1034 ; n=24) *( *24ef)(4/ pi) 
489 1035 ; = (2*Dastyet “inmghs sd/p pi) 
489 1 $ 3 = Kec, ry 4*f is an integer and C = 2*(n-24)*4/pi 
4B9 1 3 Let L = K*C modulo 8 nares =< L < 8, and let I = the integer(L) vend 
183 ! 8 : h = fract(L). then if I is even Y = h, otherwise Y = I-h 
re 1940 ; CONSTANTS: | 
459d 344 489 1906 L_INT WEIGHT = *X1E80 ; weights exponent by 61 
0489 104 “TERA _WEIG = *x1000 3; weights exponent by 32 
00000200 0489 1044 GO_TERA_WEIGHT = *x0200 ; Weights exponent by 2 in G_floatin 
00004000 0489 1045 WLMAX_WEIGHT = *x4000 3 maximum unbiasea exponent Fs F and 
0489 Hr 3 Bt dag representation 
00000019 er} + od W_ADJUST = *x19 : wg BD rat mney potnn in 
0489 1049 ; D_2_10.64: .QUAD *x0000000000006080 : 9 b4 ~ 
bcB9 1981 G_2_T0_64: 
00000000 00004410 Ree? See -QUAD %x0000000000004410 3: 2°64 
poet Wee REDUCE _LARGE : 
rea 1938 : The vores step 4 to obtain the location of the binary point in the represen- 
etd 1057 ; tation of C = 2“(n-p)*(4/pi) in two parts - the number of longwords from 
4C1 1058 ; the start and the number of bits from the most significant bit of the next 
pet) + 39 3 longword. Also K = 2*24*f must be obtained. 
52 50 *9 BF 9C 0O4C1 1963 , ROTL #-7, RO 3; Get exp field in first byte of R2 
52 19 A2 45 3 1996 SUBW He "ADJUST. 2a 3; Unbias exp and _adjust for leading 
9434 We7 ; noah R2 = location of binary 
53 9¢ 04C9 1065 ROTL  #-3, R2 R3 : pipiae’ R2 by 32 and mull by 4 to 
5 FPFFFFES ar : ise BICL #*XfFFFFFE3, R3 ; Ge of , Longwords (in bytes) + 
51 00000000'EF 4D> 68 MOVAL MTHSAL_4_OV_PI, R1 ; Get t base address of MTHSAL_4_OV_PI 
—. 3 40C 1990 SUBL R3, R1 : R} 1 points to Ist quadword of interest 
52 EO 8F tbe : 4 BICB #*XE0O, R2 : 0) = # of bits within longword 
50 7F80 8F ER 1 rg BICW #*x7F80, RO ; Clear guponent field 
50 4C€00 He 4E8 1074 BISW #*x4C00, RO ; RO = 2%2 
50 0 4ED 1075 CVTFL RO, RO ; RO=K 
us Rt | 
4FQO 1 4 : The next tep is to generate an approximation to C, call it C*" to be used 
470 1079; ” comput ng x04/ g will consist of the first three integer bits of 
4FO 1080 ; C and the first 6 beatsten bits of C. These bits will be obtained from a 
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NOTE: 


result in an inte 
so that the IV bi 
after all of the necessary fract 


MOVPSL 

BICL 

BICPSW 
Get C'* = €(0):C(1) in R3/R4. 
from the binary representation of C. 


are the the first_three bits to the 
of C(O) and bits 31:0 are the first 61 bits to the right of the binary point 


of C. 


ASHQ 
SUBL 
ASHQ 


r overflow trap. 
must be turne 


#°C<PSLSM_IV>, (SP) 


R2, (R1), =(SP) 


Be fF:28:16 


constant stored in the integer array MTH$AL_4_OV_PI. 


The ASHQ, ADDL, and MULL instructions in the follow sections may 
The overflow incurred is intentional, 
bit is not restored until 

n bits have been generated. 


; Clear integer overflow bit 


C(0) and C(1) are unsigned integers generated 
The hig 
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Generate the low 64 bits of the product 


equivalent to 
roduct is in 


At this 
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ADDL3 
BITL 
BNEQ 


More bits need to be 
not enough registers 


oint there may or m 
irst 7 fraction 
loss of significance when computin 
of significance before converting 


#*x200000, R4, 


tiplying K times C** modulo 8. 
‘ its 63:61 the octant bits, 60:24 
its of the product and 23:0 non-valid fraction bits that will be used later 
f more fraction bits need to be generated. 


. 
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enou 
s are all 1‘s or 0 
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enerated to cover the loss of significance. ; 
o hold all the potential extra bits, so that the bits 
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; Put current PSL on stack 
; (SP) = current I it 


binary point of C. 


Shift the prope 
(SP)+4 has C¢ 
oints to next quadword in 
THSAL_4_OV_PI tabl 
; Shift quadword to R3 s 
C(1) in it 


r quepeers so that 
Odin it 


e 
o that R4 has 


The result of th 


R3/R4 = K*C(1) 

Check if R4 was negative 

igned mult. error 

2 bits of K*C(0) 

L = low order bits of K*C'' 


bits in R3/R4 to 
s, there a possibi 
Consequently, we must check for loss 


If the first 8 fraction bits are 1's 
and the reduced arg 
and the reduced 


and significance will 
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3 ; already generated must be put on the stack. 


yeEey es 


tr 
>s 


SE 04 SUBL #4, SP ; Allocate an extra lLongwords on stack 
3 for later use 
6— 53 MOVa R3, (SP) ; Put L onto stack 
000006A0' EF JSB GEN_MORE_BITS ; R3/R4 contain 56 additional bits 
6— 54 ADDL R4, (SP) ; Add new high order bits to old low 
: order bits 
05 BCC ¢$ 3; Check for possible carry and adjust 
04 AE INCL (SP) ; _ high order bits accordingly 
2c BRB 4$ 7 A corry indicates a minor loss of 
: significance has occured 
14 04 AE 1D 2s: BBS #29, 4(SP), 3$ ; Check for leading ones or zeros 


; Lost significance due to leading zeros 
FFS #0, #20, 4(SP), R4 
BNEQ 
ys W°XIEFFFFFF, (SP) 
BRW LEADING_ZEROS 


3; Lost significance due to leading ones 


54 O04 AE 14 90 
6E IFFFFFFF 8F 
OOEC 


If at least one bit is set. This 
means lost significance was minor. 

One of the three high bits is set. 
lost significance was minor. 
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Octant bits are odd. Reduced argument equals 1 - f. 
BISL #*XE0000000, R4 ; Set octant bits 
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1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
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1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
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OOOO 09000909 09 09 C9 0909 09 SI SINS NN NID DD DDD DDD TTI EB ES EEE 
PAD 2 OOD NAME WIN OS WOO NAUES WIN 0 OD NAU EWN O OO NOAU EWN 0 OONOU EW O 


SOOOOOOCOOOCOOOCSOOCOOCOC COO GCOOOOOSOOOCO OOOO OOCOOCOOOOOO OOOO OOOOOOoOO 


54 04 AE 14 00 3$ FFC #0, #20, 4(SP), R4 ; If at least one bit is clear. This 
0B BNEQ : means lost significance was minor. 
6— €0000000 8F CMPL #*xE0000000, (SP) ; One of the three high bits is clear. 
0 BGEQU§ 4$ ; lost significance was minor. 
5 BRB LEADING_ONES 
53 6E 4$ MOVQ (SP), R3 ; Move bits to registers 
8E g TSTL (SP)+ ; Remove extra longword from stack 
7 
7 CONVERT: 
08 4 ADDL #8, SP 3; Clear stack 
4 : Convert bit string to double precision reduced argument. 
7 e 
4 5 Isolate octant bits and convert low order bits to double precision. 
4 EXTZV #29, #3, R4, R2 ; R2 = octant bits 
, : CVTLD 3, RO 
7 ; CVTLG 3, RO ; RO/R1 = 2°61*FLO 
? BGEQ 1$ ; Check for signed conversion error 
7 INCL R4 : _ and adjust R5 accordingly f 
: 1$ BLBC R2, 2% 3; Check for odd or even octant bits 


54 €0000000 8F 


Nate 


5 
53. 0 


53 a0 


53 

53 900 8F 
5 
03 


7E 53 
00091E80 8F 


FFEOO000 er 


Point Sine 


Ct" «te 


NMSA > > PPP PS 


DDE BEBE EEE EAI WWI IRDIPININININININIDY) 2 2 OS | QOD OOOOOOCOOOOWWOO 


~O.00 0000000000 


OOOO Cooooooooooooooo 


www taltabtabvabrabeabeah ab alah al al val eal val an val vanes 
OOPrrrrrrrrrowowownowowownwowvono 


OWT PAAOAOOO TN HANS 


MAMMA AMAIA MAUI 


& & Go oo ww ow ww www 
MOOD NAME WIN © OOD NAME WIN OOO NAME WIN MH OWO NAME WAR OO ODNAUES WN OVOONOu 


-—Oo 
reo 


SSS SSS SS SS SS SSS SSS SSS SSSSSS3Ssssseooooss 


CVTL 
ADDW 


CVTL 
ADDW 


ADDD 
SUBW 


ADDG 
SUBW 
BRW 


: Octant bits are even. 


TR TE TE TET 
> >> 
J oo 
= 
o 


© 
° 
. 
o 
. 
° 
° 
e 
oa 
° 
° 
a 
° 
e 
° 
. 
eo 
e 


LEADING_ONES: 


; bits to compute Y are available. 


4 3 
» Cosine and Sincos bas tt 74 
-SEP-1984 


R4, R3 

#U TERM_WEIGHT, R3 
R4, R3 
#GU_TERM_WEIGHT, R3 
RO, RS 

#°X9EBO, R3 


RO, R3 
#*X83D0, R3 
RESTORE 


#*XE0000000, R4 


R4, R3 
#WTERM_WEIGHT, R3 


RO, 
#*X1E80, R3 


R4, R3 
aout RM_WEIGHT, R3 
#°X03D0, R3 

RESTORE 


96:06:12 aes Macro V04-00 
29:1 MTHRTL.SRCJUVXSINCOS.MAR; 1 
3 R3/R4 = -2°298(1 = FHI) 
3 R3/R4 = “$0838 (1 - FHI) 
3 RS/R4 = =2°618(1 = FHI) 
3 R3/R4 = -2°61%(1 = FHI = FLO) 
3; R3/R4 = 1 = F 
Reduced argument equals f 
3; Clear octant bits 
3 R3/R4 = 2°29*FHI 
3 R3/R4 = 2°61*FHI 
3 R3/R4 = 2°61*(FHI + FLO) 
; R3/R4 = F = FHI + FLO 
3: R3/R4 = 2°29eFHI 
3 R3/R4 = 2°61*FHI 
3 R3/R4 = 2°61"(FHI + FLO) 
3; R3/R4 = F = FHI + FLO 


At this point it has been determined that there is a major loss of 
significance and the processing begins a looping phase. 
; the loop will generate additional extra bits of K*C' until enough significant 
During this time the stack will contain 

; four longwords allocated as follows: 


mo Shs 


Each iteration of 


SP ----> L2 : L1 and L2 are the last 2 lonquords worth of 
L1 3 signiticant bits generate 
LO : Ist three bits of LO are the octant bits 
@ : Weighting factor - keeps track of the 


number of iterations that have failed 


to produce enough significant bits. 


: If processing continues here it is known that the loss of significance is due 
; to a string of leading ones. 


MOVL 
MOVL 


LOOP_1: CMPL 


BGTRU 


R3, -(SP) 
WL_INT_WEIGHT, 12(SP) 


#*XFFEQOO00, 4(SP) 
CONVERT_1B 


3; Put L2 onto stack 
3: W = exp bias for la 
: of the product K 


; Check L2 for enough significant bits 
; Enough bits. Conve 


st_longword 
2° 


rt to floating. 


3 
VXSSINCOS ; Floating Point Sine, Cosine ond Sincos yo dt 4 1 AX/VMS Macro v04-00 Page 30 
uo08 REDUCE. LARGE r3epr 188s F881 PANANTS SResoVNStnCos.mans1 29° (38) 
QOOOO6A0'EF 16 20 : 38 JSB GEN_MORE_BITS 3 R3/R4 + cogtetag an 56 addtional 56 
P 0 
54 6€ CO 2p 1254 ADDL (SP), R4 3; Add rf order bits of existin ng 
509 1255 3 product, L2, to high order bits of 
509 1 2$ 3 _ New product, R4. 
1€ 0509 125 BCC 1$ 3 Check for corte on previous ADDL 
04 AE 06 O50B 1 28 INCL 4(SP) 3; Add carry into L1 
06 if SDE 125 BCC 1$ 3; Check for carry Mh hg begvreus ADDL 
08 A ef 1260 INCL 8(SP) ; Add carry into arry 
009 1 E3 1261 BRW CONVERT_OA ; eceeegeses to LO shmaak fraction 
5E6 1 $¢ : its are available to produce 
bee | o7 3 «equivalent argument 
04 AE FFFFFFFF 8&F D1 bzee 1265 1$: CMPL #-1, 4(SP) 3 Check for L1 all 1's 
1c =A ef : 66 BGTRU CONVERT_1A 3 tee 1's es Enough precision bits 
; © compute Y 
6— 53 7D ber 1268 mMOvVa R3, mpress representation of L 
FFC8 OC AE 1000 8F 4000 BF 3D OSF3 1 $2 ACBW #_ max PMEIGHT, #W_TERM_ uéignre 12(SP), LOOP_1 
OSFE 1270 : Increment weight tag factor. If 
OSFE 1271 3 wg iahting factor is greater than 
OSFE 1 Le: 3 128 then no more bits need to be 
OSFE 127 3 = generated. 
Oere 1978 
OSFE 1276 ; The es factor is greeter than 128. This means that the reduced 
OSFE 1277 ; eroyeen is either not q [ot toguighesle from 1 or too small to be represented 
OSFE 1 4 3; in F_format (i.e. underflow.). Zero is returned in R3 for the reduced 
OSFE 1279 ; argument to signal this occurance. Note that under these conditions the 
OSFE 1280 ; correct function value is one of the values 0, +/-1. The 
OSFE 1281 ; correct choice is determined by the calling program based on the octant bits 
OSFE 1 Hf 3 returned in Ri. 
OSFE 1283 ; 
53 D4 OSFE 1284 CLRL R3 ; Reduced argument is zero 
52 08 AE O03 1D EF Q600 1285 EXTZV #29. a3. 8(SP), R2 3; R2 = octant bits 
5E 10 CO 0606 1286 ADDL 3; Clear stack 
008E 31 0609 1287 BRW RESTORE 
bebe 1288 
be C 1289 CONVERT_1A: 
6—E 53 DO 060C 1290 L R3, (SP) 3; Put L2 onto stack 
O60F 1291 CONVERT _1B: 
Bene 1 38 5 
red : 32 3 CVTLD (SP)+, RO 3; RO/R1 = 2*W*L2 
50 BE “erp O60r 1295 CVTLG (SP)+, RO 3; RO/R1 = 2*W*L2 
1 61 1 38 BGEQ 1$ ; Check for signed conversion error 
6— 06 $2 ' oe INCL (SP) 3 and adjust L1 accordingly. 
° : 
0817 1299 ; : 1$: CVTLD (SP),R3 
BI? 330 | 
53 6€ 4EFD $e15 1 ¢ 1$: CVTILG (SP), R3 ; Convert L1 to double and check for 
07 ~=«18 718 ' ? BGEQ 2s 3 signed conversion error. 
ef ' i : ADDW #*x9000, R3 
53 8200 AQ 061D 1307 ° ADDW #*x8200, R3 : R3/R4 = 2%W*(2464 = L1) 
Bt 9 $22 1308 BRB 3$ 


B 4 
VXSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 02:06:1 AX/VMS Macro V04-00 Page 31 
Oe REDUCE. LARGE 66-SEP-1984 96:98:18 EMTHRTL. SREJOVXSINCOS MAR; 1 . 30) 
$3t 1840 
624 1311 ; 28: ADDW #1900 R3 
624 1 1g ; SUBD3 R3, 6_2.10_64, R3 
e54 j314 * 
53. 0200 8F AO 0624 1315 28: ADDW oO 208. R3 3 R3/R4 = gaust] 
53 FEBA CF 53 43FD be 4 ! 1 SUBG3 R3, G_2_T0_64, R3 3 R3/R4 = 2*%W*( 2°64 = L1) 
6 1318 ; 
6 5 1 i§ ; 3$: SUBD RO, a 3 R3/R4 = 2*W2(2°64 = L1 = L2) 
6 8 1320 ; EXTZV #29, #3, (SP)+, R2 : Re = octant bits 
630 1321 : SUBL SPS+, R3 > R3/R4 = Y 
06 0 1 : ; BRW RESTORE 
3 0 : 4° 
53 50 42FD 3 0 1325 3$ SUBG2 RO, R3 3 R3/R4 = 2*W#(2%64 = L1 = L2) 
52. BE 03 if . EF 06 4 6 EXTZV #29. #3, (SP), R2 : R2 = octant bits 
gg Hi ria $63 1328 SUBG2  (SP)+, R3 : R3/R4 = 
0056 31 0641 1 $3 BRW RESTORE 
get 1 
0644 1 3 LEADING_ZEROS: 
Sit ee 
0644 1335 ; If processing continues here it is known that the loss of significance is due 
0644 1336 ; to a string of leading zeros. Note that it is known that the loop for 
0644 1337 ; leading zeros will terminate before an underflow condition occurs so that the 
Bete 1336 3; loop does not include a test for underflow. 
7E 53 DO 0644 1329 MOVL R3, -(SP) : Put L2 onto stack 
OC AE OO001E8O BF DO 0647 1341 MOVL #L_LINT_WEIGHT, 12(SP) ; W = exp bias for Last longword 
Rear ! 4g : of the product K*C' 
04 AE OOTFFFFF 8F 01 O64F 1344 LOOP_O: CMPL #*XOO1FFFFF, 4(SP) ; Check L2 for enough fraction bits 
19 0657 1345 BLSS CONVERT_O : Enough bits. Convert to floatin 
OOOOD06A0'EF 16 pees 1347 JSB GEN_MORE_BITS 3 sae sayrem an addtional 56 bits 
a © ' 
54 6€ CO Q65F 1348 ADDL (SP), R4 ; Add low order bits of existing | 
0662 1349 ; product, (SP), to high order bits 
0662 1350 3 _ of new product, R4. 
03 1€ 0662 1351 BCC 1$ ; Check for carry on previous ADDL 
04 AE 06 b66s 1 26 INCL 4(SP) ; Add in carry. Note: High bits of R4 
67 135 ; are 0 so no gerry possible on INCL 
54 DS 0667 1354 1$ TSTL R4 3; Check for all O's ; 
12 0669 1355 BNEQ CONVERT_OA ; Not all O's. Enough precision bits. 
6E 7D 0668 1 2$ MOVQ ° ; Compress representation of 
OC AE 1000 8F AO $93 135 ADDW #W_TERM_WEIGHT, 12(SP) ; Increment weighting factor. 
09 =i or : 38 BRB LooP_0 
676 1360 CONVERT_OA: 
6— 54 00 078 1361 MOV R4, (SP) 3; Put L2 onto stack 
679 136¢ CONVERT_OB: | 
679 1364 : CVILD (SP), RO ; RO/R1 = 2*weL2 | 
679 1365; 
| 
| 
| 


. 
C 
mg : Floating Point Sine, Cosine and steeee aes b Oe: 96:12 ale Macro V04-00 Page 36 | 
-00 REDUCE_LARGE =SEP-1984 UMTHRTL. SREJUVXSINCOS..MAR; 1 ( i 
50 g& = 679 1 66 CVTLG {3P). RO 3 RO/R1 = 2*We*L2 
Oe "18 067 136 BGEG i If L2 is negative, adjust L1 to | 
6— 06 7 oe rt: INCL (SP) ; to reflect signed rather than 
? ! ! $3 3 unsigned conversion | 
681 1371; 18: CVTLD (SP) ; Convert L1 to double | 
681 1 ie : ADDW ow lemme: “WEIGHT, R3 : Rg = O*WeLt 
681 13735; ADDD RO; : R3/R4 = 2*We(L1 + L2) 
681 1374 : EXTZV #29, #3,_(SP)+, R2 : Re = octant bits 
6 0 : SUBL  (SPS+, re : = Y 
5 ge 4EFD 0681 1 7 1$ CVTL (SP) 3; Convert L1 to double 
53 9 00 AO 0685 1 28 ADDW2 #Gw_f Re _WEIGHT, R2 ; Rg7R = *we 
50 40FD 68A 137 ADDG RO, R : R3/R4 = 2*We(L1 + L2) 
52 8 O03 1D EF O68E 1380 EXTZV #29, #3, (SP)+, R2 : R2 = octant bits 
3 BE GEFD 0693 1 81 CVTILG (SPS+, =(SP) 
5 8E C2 0697 1 ¢ SUBL2 (SP)+, 3; R3/R4 = 
Seon 1384 | 
069A 1385 RESTORE: 
6— B8 O69A 1 36 BISPSW (SP) ; Restore IV bit and exit 
SE 04 CO O069C 138 ADDL #4, SP : Remove mask from stack 
05 +444 1338 RSB 
6A0 1390 GEN_MORE_BITS: 
Sead 1392 
06A0 1 38 : This subroutine generates 56 extra fraction bits and stores them in R3/R4 
SS 1s 
04 ce 06A0 1396 SUBL #4, R1 ; Adjust pointer to next quadword 
53 61 52 79 O6A3 1397 ASHQ R2. FRI). R3 : Put next coefficient in R4 
53 6000 50 7A Q6A7 1398 EMUL RO, R4, #0, R3 ; Generate the fraction bits in R3/R4 
03 18 O6AC 1399 BGEQ 1$ : Check for signed arithmentic error 
54 50 CO O6AE 1400 ADDL RO, R4 ; Negative factor, adjust R4 
05 0681 1401 1$: RSB 
pees 1406 
06B2 140 
06B2 1404 
06B2 1405 


nt Sine, Cosine and tinees i ab tae 9 Oe: 96: 1 ies 8 o V04-00 Page PY | 


mtg : looting Poi 
-00 REDUCE _DEGREES EP-1984 MTHRTL.SRCJUVXSINCOS..MAR; 1 1) 
p68 1407 -SBTTL REDUCE_DEGREES 
68 1209 : This routine assumes that the absolute value of the 7proyaent is in RO. 
6B2 1410 ; The reduction process is performed in two grages. The first stage of 
6B2 1411 ; the reduction reduces the argument modulo 360 to a value P as that 2°24, 
re 1ei¢ ; and the second stage reduces the argument to less than 4 
06B2 1414 ; Constants used in this reduction: 
pe8 1415 ; 
68 1218 
06B2 1417 POWER_MOD 360_0: ; Powers of 2 modulo 360 for t1 = 0 
0008 0004 4 3 poor 68 1213 -QORD 1» ° 4, 
0080 0040 00 010 O6BA 141 «WORD 16, 32, 64, 12 
OOF8 0130 0098 0100 pete 12 y WORD 256, 152, 304, 248 
O6CA 14 ¢ POWER_MOD_360_1: eee, gf 2 modulo 360 for ti <> 0 
0008 00B8 0110 0088 O6CA 14 QORD” 136, 572 184, 
0080 0040 0020 0010 0602 1424 WORD 16, 64, 128 
OOF8 0130 0098 0100 pepe 1 4 WORD 256, 1 304, 248 
O6E2 1427 REDUCE_DEGREES: 
50 4C€80 8F B1 Q6E2 1428 CMPW #*x4C80, RO : Compare ix: with 2°24 
14 Q6E7 1429 BGTR LAST T_STEP ; Branch to special logic for med arg 
06E9 1430 
06E9 1431 ; 
BOE® 1336 3; It is assumed here that the argument is greater than 2°24. 
06E9 1232 ; The argument is reduced as follows: 
06E9 1435 ; x “tef, where t > 24 and 1/2 =< f <1. And let J = 2*°24*f = 
06E9 1436 ; 3815eu1 + J2 > K = = §4 (t- a Since 2°15 = 8 modulo 360, we have that 
06E9 1437 ; J = 8*J1 + J2 modulo 360. ‘st - 24s in. + t2. Note that 
06E9 1438 ; (2*12)*2 = 2°24 = (Ba9)a(281 ) = (2*9)#(2*3) = 2°1 - dulo 369. _Hence, 
06E9 1439 ; if tl_is_ not_zero, K = * = 2*°(12eti+t2) = C2*(12* t1 )J* os 
06E9 1440 ; C2*123*C2*t2) = 1$682*t agéy\o 36 For ti = 0, K = alle 
O6E9 1441 ; define K' congryent to 2*t2 if t1=0 and gongrugns to 1 6294t2 otherwise, 
06E9 1346 3 where O=< K' Then x" = K'*(8*J1 + J2) is congruent to x modulo 
06E9 1443 ; 360 and x" < 2* 
06E9 1444 
943 1445 
51 50 _Q00007F80 8F CB O6E9 1446 BICL3 #*X7F80, RO, R1 H 
51 4C00 8F A8 O6F1 1447 BISW #*x4C00, Ri ; RI = J 
51 C2 bere 1223 SUBL R1, £0 ; RO = t'*2*7 
52 51 7FFFOOOO He 8 pery 1450 BICL3 #*X7FFFOOOO, R11, R2 : R2 = J1*2°15 
51 2 4 01 1451 SUBF R2, R 3; RI = J2 
52 0600 8F ac 4 1e06 SUBW Pee R2 ; R2 = J1 
FS ae 4 , 1¢37 ADDF R2, ; R1 = 8J1 + J2 = J modulo 360 
50 50 F9 8F 9 soe 1455 ROTL #-7, RO, RO ; RO=t' 
52 50 OC 7 0711 1596 DIVB3 ig. Re R2 ; R2 = tl 
52 96 4 0715 145 MULB ; Re = es 
50 2 Bo i8 1oe8 SUBB R2, ; RO=t 
36 95 0718 1460 TSTB R2 ; Check for tl = 0 and choose Kk’ 
0 12 071D 1461 BNEQ 1$ 3; _ accordingly 
50 8F AF4O 4D O71F 1468 CVTWF POWER_MOD_360_OCROJ, RO ; RO = kK’ 
05 11 724 146 BRB 2$ 


Er Res 


E 4 
VXSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-19 71 AX/VMS v04-00 Page 34 | 
ByaO8 REDUCE. DEER aes 94: 98: 18 YATARTL. SUvX SINCOS.MAR; 1 t) 
50 a0 AF4O 4D 726 4 1$: CVTWF POWER MOD_360_1CROJ, RO ; RO = 
5 51 44 4 3 5 2$: MULF R1, RO 3; RO = x (mod 45) 0 =< RO < 2°24 
tSE 1269 
4 3 8 LAST_STEP: 
f : y ; Argument reduction scheme for arguments with absolute value less than 2°24 
oF E gi : The reduced argument Y is computed as follows: 
72E 5 Let 1 = 
072E 4; I is even 
Bf E S$: then Y = X = $371 
4 : § 5 else Y = (1#1)#45 = x 
53 50 B86 8F 0B603DB6 gf 54 oF 3 8 EMODF #LF_1_0V_45, #X_1_0V_45, RO, R3, R4 
0739 9 ; R3 = I = integer part of ix:/45 
1653 €9 0739 0 BLBC RS. EVEN 3 Branch if octant bits are even 
54 53 (01 C1 073C 1 ADOL3 » R3, RS 3; RO 2141 
54 54 4E 15 $ CVTLF 4 3 
54 00004334 BF 44 074 MULF oo R4 > RG = 45*(14+1) 
54 50 42 O74A 4 SUBF . 7 RO = Y 
53 F8 8F 8A 074D 5 BICB #*xXF8, R3 ; Save only Last three octant bits 
05 oral § RSB 
54 53. 4E 0732 8 EVEN: CVTLF R3, R4 ; RO =] 
54 0000C334 8F 44 075 ) MULF #LF_M45, R4 3 RG = -45e] 
54 50 40 075C 0 ADDF . RG = Y 
53 F8 BF 8A 075F 1 BICB #*xXF8, R3 3; Save 
05 0763 ¢ RSB 
0764 


| 
only last three octant bits 


DOOOWOOOOOCOCOOOOO W-- 


1 + jt 


F 4 
nt Sine, Cosine and Sincos 16-SEP-1984 1 AX/VMS v04-00 P 5 
aes 4: 96: 18 YATARTE og SUvX SINCOS.MAR;1 ” (33) 


»SBTTL RADIAN_POLYNOMIALS ; Polynomials for arguments in radians 


Polynomial evaluation for COS(Y) for Y in radians 


wo Sete te 


COS_R: 


OWONAUE WN OOONOM 


; Polynomial evaluation for -COS(Y) 


N_COS_R: 


ing P 
DEGREE 
764 14 
764 14 
764 14 | 
764 14 
it | 
ree 48 | 
764 15 
764 1 
764 1 _COS_R: 
53 8000 8F AA 0764 1 eicy #*x8000, R3 ; R3/R4 = {Y! 
769 «#1 3 CMP #*x4000, R3 Compare 1/2 with {Y! 
4 : 3 Using CHPL for G does not yoeld enough signif fraction bits 
50 53 33FD $769 12 CVTGF R3, RO ; Need to conform to old code. 
50 4000 8F B61 8768 1510 CMPW #*x4000, RO : Compare 1/2 with iY} 
2—E «618 O772 1511 BGEQ LEQL_HALF ; Sufficent overhang is available 
0774 131g NEEDS_DOUBLE : 
0774 1513; 
0774 1514; MULD R3, R3 3 R3/R4 = YeY 
0774 1515; 
0774 1516 
53 53 44FD 0774 1517 MULG2 R3, R3 3 R3/R4 = YeY 
0778 1518 
0778 1519 NEEDS DOUBLE SINCOS: 
0778 1520 ; PUSHE ; Save high half of Yev 
7E 53 7D 0778 1521 mMOva R3, -(SP) ; Save G rep of Y* 
54 FFF 8F AA O77B 15 ¢ BICW2 #AXIFFF, R4 ; Zero the rounding bit (bit 12) and 
0780 15 ; _all lower order bits of word 4. | 
53 53 33FD th 44 13 : CVTGF §R3, R3 3; R3 = Y*2 truncated. | 
FRAG CF 04 53 55 0784 1526 POL YF R3, poet 7 COSTBR2; RO = An} hi 
078A 1527 ; SUBL3 #*X80 +, RS > R3/R4 = YeY/ 
078A 1528 ; SUBD = #1, RS > R3/R4 = y275 = 1 
078A 1368 : SUBD R3, RO 3; RO/R1 = COS(Y) 
078A 1530; CVTDOF RO, RO ; RO = COS(Y) 
4 1531 
53 8€ 10 C3 O78A 1336 SUBL3 #*x0010, (SP)+, R3 3 
54 «BE DO O78— 15 MOVL (SP)+, R4 3 R3/R4 = Y*2 
53 08 42FD «0791 «+1534 SUBG2 #1, R$ : R3/R4 = Y*2/2 = 1 
50 29 99FD 0795 139? CVTFG RO, RO 
50.) (65 a4 799 1536 SUBG2 R3, RO 3; RO/R1 = COS(Y) 
50 50 33FD O079D 1537 CVTGF RO, RO 3; RO = COS(Y) 
05 O7A1 1 } RSB 
TAZ 1539 LEQL_HALF: 
7A2 1540 ; CVTDF R3, R4 3 RG = Yt 
54 53 33FD O7A2 1541 CVIGF R3, R4 7 RG = iYt 
54 54 4h 1366 MULF R4, RS 3; R4 = YY 
F871 CF 03 54 £55 7 154 POL YF R4, #COSLENR1-1, COSTBR1; RO = COS(Y) 
05 1544 RSB 
136 
7 1 rt 
154 
1268 
154 
1550 
1551 | 
| 


-=cC 


—— — ———, 


| 


6 4 
XSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 06:1 AX/VMS Macro V04-00 Pa 
Wao RADIAN_POLYNOMIALS ; Polynomials for arg ats P 96:98:18 EMTHRTL SREIOVXSINCOS.MARS 1 - 3%) 


53 8000 8F AA 0780 1 § BICW 9. 18000. R : R3B/R4 = fY} 
785 (1 : CMPW #*x4000, R ; Compare 1/2 with Y: 
i i ¢ : Using CMPL for G does not yoeld enough signif fraction bits | 
50 53 33FD 0785 1 ; CVIGF §=6R3, R ; Need to conform to old code. 
50 4000 §F 8) 9789 1 5 CNPue Re Ka000, RO S Compare 1/2 with 1Y; | 
—E 18 O7BE 1 BGEQ 1 ; pySiicent overhang is available 
7C¢ 1559 ; MULD R3, R3 ; RS/R4 = YY 
4s : $9 : PUSHL R ; Save high half of Y*Y 
rE ; 44FD 7¢ 1368 MULG2 Rie R3 3 R3/R4 = YeY 
E 7D rts 130 MOVQ R3, -(SP) 3; Save G rep of Y*2 
54 FFF BF AA OQ7C7? 1564 BICW2 #*XIFFF, R4 ; Zero the rounding bit (bit 12) and 
7CC 1392 3 Ll lower order ye of word 4. 
53 53 33FD re 1366 CVTGF R3, R3 ; R3 = truncated Y* 
FBSA CF 04 53 55 708 1308 POLYF R3, #COSLENR2-1, COSTBR2; RO = Y*4*Q(Y*rY) 
bape 1390 ; SUBL3 #*x80, (SP)+, R3 3s R3/R4 = Y¥eY¥/2 
706 1571; SUBD ss #1, RS : R3/R4 = -(1 = ¥*2/2) 
07D6 1376 : SUBD RO. R3 : R3/R4 = -COS(Y) 
Bope 1573 ; CVTDF R3, RO ; RO = -COS(Y) 
Bee 1s 
53 8E 10 C3 0706 1378 SUBL3 #*x0010, (SP)+, R3 3 
54 8E DO O7DA 157 MOVL (SP)+, R4 3 R3/R4 = Y*2/2 
35 90 Sore Oven 1879 pase #). © eae peodke- 
28 50 Ger O7E5 1580 SUBG2 RO, R3 : R3/R4 = -COS(Y) 
50 53 33FD O7E9 1581 CVIGF = R3, RO ; RO = -COS(Y) 
O7ED 1386 
05 O7ED 158 RSB 
OrEE 1888 
Ore 1336 ; 1$: CVTOF R3, R4 3 RO = {YI 
tae a re: 
F825 CF 4003) 54S) sOO7FS «1590 POLYF R4. #COSLENRI-1, COSTBR1: RO = COS(Y) 
50 8000 8F AC Orr 1591 XORW #*x8000, RO ; RO = -COS(Y) 
05 0800 1336 RSB 
O80 1894; 
0801 1399 : Polynomial evaluation for -SIN(Y) 
» % 36 : | 
BT a sie 
53 8000 8F AC $801 1399 ~~" KORW «= #*X8000, R3 
BS 1 
a 1608 ; Polynomial evaluation for SIN(Y) 
1604 ° 
Og 1605 P_SIN_R: 
06 16 $ ; PUSHL R : Save high half of Y 
7E 53 7D 38 198 MOvVaQ R35, -(SP) ; Save R37R4 = Y- 2 E 
54 FFF 8F AA 1608 BICW2 #*XIFFF, R4 ; Zero the rounding bit (bit 12) and 
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Olynomials for arg 6-SEP-1984 
Re R3 
RS gp BSINGENR1, SINTBR 
ie 
RO, R 
(S$P)+, RO 


R 
(phos, 80 


H 4 04: 
feet fee esnees erie lake Teele Daten. SREJUVXSINCOS.MAR: 1 © 


git lo 
R5 = : 
Save Y 


= Ye 


= P( 


RO = tr 
ee 


X/VMS Macro V04-00 


wer order bits in word 4. 
truncated. 
eet ot 


Y 
Y*y) 


uncated Y * P(Y*2) 
a 


m 


-C 


I 4 
XSSINCOS ating Point Sine, Cosine and Sincos 16-SEP-1984 02:06:1 AX/VMS Macro v04-00 Page 38 
yyo08 * POLYNOMIALS ; Polynomials for argu ~$eE 138 99398318 UMTHRTL « SREJUVXSINCOS MAR: 1 . (38) 
4 16 8 -SBTTL CYCLE_POLYNOMIALS ; Polynomials for arguments in cycles 
D 1631 
gia 
p 18 4 ; Polynomial evaluation for COS(Y) for Y in cycles 
D 1636 - 
D 16 $ P_COS_C: 
D 16 8 3 CMPF #LF_2. OV_PI, R3 3; Compare 3/pi with {¥! 
53 00000000 5F304004 8F ds D 16 CMPG #G_2_Ov_PI, R3 3; Compare 2/pi with iY: : 
a4. 9 1969 BGEQ LEGL-2_OVv_PI 3 syiticent overhang is available 
1641 ; MULD R 3 RS/R4 = YeY 
8 1966 3 PUSHL 3 Save high half of Yey 
53. 53 44FD 1644 MULG2 R3, 3; R3/R4 = YeY 
7E 5 7D O83F 1645 MOVQ R35, -(SP) 3 SAVE Y*2 
54 FFF 8F ots 1946 BICW2 #*XIFFF, R4 ; Zero the rounding bit (bit 12) and 
oat 164 : all lower order bits of word 4. 
53 33 33FD vt 1648 CVTGF R3, R3 ; R3 = Y*2 truncated. 
F817 CF 04 55 0848 1649 POLYF R3, #COSLENC2-1, COSTBC2; Rn = ¥*42Q(YrY) 
0851 1650 ; SUBL3 Mele: & (SP)+, R3 3 R3/R4 = Y¥RY/4 
0851 1651 ; UBD #1, R 3 R3/R4 = Y*2/4 = 1 
+3 163¢ 3 UBD R3, RO 3; RO/R1 = COSC(Y) 
+3 1937 3 CVTDF RO, RO 3; RO = COS(Y) 
53 8E go C3 0851 1655 SUBL3 #*x0020, (SP)+, R3 3 
54 f DO 0855 1936 MOVL (SP)+, R4 3 R3/R4 = Y*%2/46 
age Se ee Feseticak 
50.53 ser 0860 1889 SUBG2 R3, RO : RO/R1 = COS(Y) 
50 503 A} Fe) 1969 eyiGr RO, RO ; RO = COS(Y) 
bpeo 1956 LEQL_2_0V_PI: 
0869 1663 ; COUTOF R3, R4 3; RO Vt 
a Ee bbe tees MULE, Ra! Re te = ve¥ 
F7DE CF 0K SKC 0890 1666 POLYF R4. MCOSLENC1-1, COSTBC1: RO = COS(Y) 
05 pe7g 1667 RSB 
0877 1668 
Bre 
0877 1671: Polynomial evaluation for -COS(Y) 
77 1676 : 
NER cos. 
a4 1675 :~ ~ CMPF #LF_2 OV_PI, R3 : Compare 2/pi with {¥} 
53 00000000 5F304004 8F SIFD 77 1676 CMPG #G_2_6v_PI, R3 3; Compare 2/pi with ‘Yi Z 
ee «C18 Ht 167 BGEQ $ : Sufficent overhang is available 
1978 MULD Re R3 3 R3/R4 = YeY 
1927 ; PUSHL R ; Save high half of Y*yY 
53 33 44FD 1681 MULG2 RR R3 : R3/R4 = Y*2 
7E 3 7D 9 16 ¢ MOVa R35, -(SP) 3; Save Y*2 
54 FFF 8F AA £ BICW2 #*XIFFF, R4 ; Zero the rounding bit (bit 12) and 
1 1686 : gt lower order pits of word 4. 
53 53 33FD 91 1685 CVTGF = 6R3, R3 ; R3 = truncated Y“2. 
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J 4 
nt Sine, Cosine and Sincos 16-SEP-1984 02:06:1 AX/VMS Ma v04-00 Pp 9 f 
ALS ; Polynomials for argu aes 99:98:18 ERTMRTL. SRETOVXSINCOS.MAR: 1 a Bay | + 


v 
°o 


POLYF &3, #COSLENC2=1, COSTBC2; RO = Y*4*Q(Yey) 
; SUBL3. #*k100, (SP)+, R3 + R3/R4 = YR¥/4 
: SUBD #1, R : R3/R4 = -(1 = ¥*2/46) 
: SUBD RO, R + R3/R4 = -COS(Y) 
; CVTDF =6ORS, RR : RO = -COS(Y) 
SUBL3 #*x0020, (SP)+, R3 ; 
MOVL (SP)+, Ra + R3/R4 = Y*2/4 
SUBG2 #1, RS + RB/R4 = <(1 = Y*2/4) 


CVTF R0. RO 


SUBG RO. R ; R3/R4 = -COS(Y) 


CVIGF RO : RO = -COS(Y) 
RSB 

- 1$ CVTDF 3. RG : RO = ft! 

is: CVIGF R3, R4 + RG = {YY} 
MULF 4, R4 > RG = YeY 
POLYF R4, #COSLENC1-1, COSTBC1: RO = COS(Y) 
xORWe #*x8000, RO : RO = -COS(Y) 


; Polynomial evaluation for -SIN(Y) 


SINC: 
XORW #*x8000, R3 ; R3/R4 = - Y 


: Polynomial evaluation for SIN(Y) 


PI tt nd dd ed ed ad ed md DOODODOOOOCOOOOOOOOOOWWOc 


RO ODNOA NE WN = OOO NAUE WIN (OOO NA UE WIN 9 OD NAU SW 0 OD NOU EWI ODN 


P_SIN_C: 
MOVQ R3, -(SP) 3; Save argument 

: MULF R3, R3 
BICW2 #*XIFFF, R4 ; Zero the rounding bit (bit 12) and 

3 alt lower order bits of word 4. 

CVTGF R3, R3 ; R3 = Truncated Y 
PUSHL R53 ; Save Y truncated 
MULF 2 Ree R3 3; RB = Ye 
POLYF R #SINLENC-1, SINTBC ; RO = PC(Y*Y) 

: MULF_ (SP), RO : RO = YeP(yey) 

: SUBL3 #*x100, (SP), R3 : 

F} MOVL 4(SP), R4 3s R3/R4 = Y/4 

: SUBD rs R3 > R3/R4 = 3/42Y 

; SUBD =—saR’3,,_-RO + RO/R1 = SIN(Y) 

; CVTOF RO, RO ; RO = SIN(Y) 
MULF (SP)+, RO ; RO = truncated Y * P(Y*2) 
SUBL #*x0020, (SP), R3 : Form the quent ity 
MOVL 4(SP), R4 > R3/R4 =_Y/4 
SUBG2 (SP)+, R3 : R3/R4 = 3/42Y 
CVTFG RO, RO 

74 SUBG R35, R 3 RO/R1 = SIN(Y) 


3 aa ivELBMBSOMGRIRES® Mba fess eGos ness TECREECHBES 99:06:18 WAKE. SRESEARta@Be man;r POPP Rs 


50.0550 3560 + | $3 gyicr RO, RO ; RO = SIN(Y) 
Fe 1428 er 


L 4 
VXSSINCOS Floating Point Sine, Cosine and Sincos 16-SEP-1984 :06:1 AX/VMS Macro V04-00 Page 41 
ett BEGREE. POLYNOMIALS pate 7 94:98:18 ERTHRTL. SREJUVXSINCOS.MARS 1 9° (35) 
FE 1747 .SBTTL DEGREE_POLYNOMIALS | 
nS | 
FE 1750 P_COS_D: 
54 2EET4ZES BF $1 OBEE 175] ~~ CHPE ALF_90_OV.PI, RG ; Compare 90/pi with Y | 
a: 8 17 § BGEQ 2 3 Double precision isn't needed 
rime: a OB RS he 
3 u 17 5; PUSHL Ra’ : Save high half yey 
53 4 99FD 0907 17 5 CVTFG R4, R : R3/R4 = Y 
E : nas 3 F 5 § move 4 * =(SP) ; rahe y*2 
54 so 1FFF BF gi¢ 1760 BICW2 #*XIFFF, R4 : Zero the rounding bit (bit 12) and 
91 1761 5 Ll lower order bits of word 4. 
53 53 33FD a7 1766 CVIGF R3, R3 : = Y*2 truncated. 
F778 CF OS BB ook ont 1988 we ee ee ee bith ie 
0925 1766 ; SUBL3 #*X680, (SP)+, R3 : RS/RG = ye¥/2"13 7 
09 3 1767 : CMPL © #*X1000, R4& : Check for loss of significance when 
53 BE 00000000 BF c3 0925 1769 SUBL3 #*x00D0, (SP)+, R3 ; 
54 BE 00 092d 1770 MOVL +, RG ; R3/R4 = Y*2/2°13 ! 
54 00000200 8F 01 bazo 44 CMPL #*x0260, R4 : wpeee See weet of stent ticance when 
0937 1558 : D_ float used *x1000. This converts to 
0337 HAs : *x0200 in G_floating. 
12 1A 0937 1776 BGTRU 1$ ; 
7E 53 7D 0939 1777 MOVQ = R3,_ = (SP) : (SP) = ¥*2/213 meek 
at a a OR 
093C 1780 : ADDD = (SP) +, RO : RO/R} = r(Y¥*2) + Low bits of ¥*2/2°13 
093C 1781: 1$: SUBD #1, R$ + R3/R4 = -(1 = YeY/2*13) 
093¢ 1788 ; SUBD = R3,_ RO : RO/R1 = COS(Y) 
8 : 1784 3 CVTDF RO, RO ; RO = COS(Y) 
54 1FFFOOOO 8F CA Soa 1785 BICL2 #*X1FFFOOOO, R4 ; R3/R4 = high 40 bits of ¥*2/2°13 
6— 53 42FD 0943 1786 SUBG2 3, (SP) ; (SP) = low 13 bits of Y*2/2*13) 
50 BE 4OFD 0947 1787 ADDG2 (5P)+, RO ; RO/RI = r(y*2) + low 13 bits of ¥*2/2*13 
5 08 42FD 0948 1788 1$: SUBG2 «#1, RS > R3/R4 = <(1 = ¥*2/2°13) 
50 3 42FD O094F 1789 SUBG2 R3, RO : RO/R1 = COS(Y) 
50 50 33FD 09 : 1790 CVTIGF RO, RO > RO = COS(Y) 
957 179 
05 997 1798 RSB 
54 54 hh 338 1794 2$: MULF2 R4, R4 > RG = YeY 
F728 CF «6003854 SCS (0958 «1795 POLYF R4, #COSDLN1, COSDTB1 ; RO = COS(Y) 
05 0961 1796 
as 8h 
962 1799 N_COS_D: 
54 2€€142€5 ef 51 0962 1800 ~ ~ CMPF  #LF_90_OV_PI, R4 ; Compare 90/pi with Y. 
1 18 0969 1801 BGEQ ; Double precision isn't needed 
968 1 Oe ; CVTFD R4, R ; RS7RG = Y 
968 1 : MULD R3, R 5 = Yey 


M 4 | 
XSSINCOS ; Floating Point Sine, Cosine and Sincos 16-SEP-1984 71 AX/VMS Ma v04-00 FP 
pet BEGREE. POLYNOMIALS patie oF: 96: 8 YRTRTL. SRCJUVXSINCOS MAR: 1 ek 3 : 
368 } : : PUSHL R3 ; Save high half of yey | 
4 99FD Q9 1 CVTF R4, R 3 R3/R4 
33 3 46rd 96F § HULG2 RS, RR i REIRG = Y*2 
E 7D 0973 1 8 MOVa R3, -(SP) : Save R3/R4 = Y*2 
54 FFF 8F AA are | o? BICW2 #*XIFFF, R4 : core pee apace bn A 44 cote i, and 
P r i 
53 5333D 0978 1811 CVIGF 3, RB ae oe ig tamed 
F717 CF 9 33 5 97F 1 ig POLYF R3, #COSDLN2, COSDTB2 3 RO = r(¥e 
0 0 99FD 0985 1814 CVTF RO, R : Need RO in v$9 ft ating format. 
989 1815 ; suaL3  #°K680, (SP)+, RS : RSIRG = 
2 4 ' 1$ ; CMPL #*x1006, R4 : Check Ry is,’ of significance when 
53 BE 0000900 Yl C3 0989 1 18 SUBL3 #*x00D0, (SP)+, R3 $ 
4 8€ DO 0991 181 MOVL (SP)+, R4 t RB/R4 = ¥*2/2°13 
54 00000200 8F D1 0994 1820 CMPL #*x0260, R4 3; Check for loss of stent ficance when 
998 1821 3 subtracting 1 1358." us 
998 1 § 3 dD. float ysed *x1000. T Ys converts to 
Bees } 7 : ™x0200 G_floating. 
12 1A 0998 : 2 BGTRU 1$ :. 
7E 53 7D 0990 is 7 MOVa R3, -(SP) : (SP) = ¥*2/2°13 
QO9A0 18 3 5 BICL #*XFFFFOOOO, R4 ; R3/R4 = high bits of UR ye mh 
O9A0 \3 j SUBD R3, (SP) ; (SP) = Low bits of Y 
O9A0 1830 ; ADDD (SP)+, RO 3: RO/R1 = r(¥*2) + {Sy bits of ¥*2/2*13 
09A0 1831 > 16 SUBD #i, RS t R3/R4 = <(1 = Ye¥/2°13) 
O9AD 1 3 3 SUBD RO, R 3; R3/R4 = -COS 
Hy 1834 3 CVTDOF R3, RO 3; RO = -COS(Y) 
54 1FFFOOOO 8F CA OQO9A0 1835 BICL #*X1FFFOOOO, R4 3; R3/R4 = high 40 bit f ¥*2/2*°13 
6E 53 42FD O9A7 18 § SUBG R3, (SP) 3 hye = lo W713 bits ot y" *2/2°13 
29 8E 40FD 9AB 1 ADDG (SP)+, RO ; /R1 = rty*2) + low 13 bits of ¥*2/2*13 
53 08 42FD ~O9AF 1838 1S: SUBG2 #1, RS : R ny = =(1 = ¥*2/2*13) 
53 50 42FD 0983 1839 SUBG RO, R3 3; R3/R4 = -COS(Y) 
50 53 33FD 0987 1840 CVIGF = R3, RO : RO = -COS(Y) 
98B 1841 
05 4 ' ts RS3 
54 54 44 apes 1844 2$: MULF2 R4, RG 3; R4 = YeY 
F6C7 CF 03 54 55 SBF 1845 POL YF R4, #COSDLN1, COSDTB1 3; RO = COS(Y) 
50 8000 8F AC 9¢c5 1 #8 XORW #*x8000, RO 3; RO = -COS(Y) 
05 SCA 184 RSB 
ae is 
9CB 1850 N_SIN_D: 
54 54 52 ace ; 23 > sin 8 MNEGF R4, R4 3 R4 = -¥ 
50 654 + i3 O65 ' 38 pet MULES _ _ RO 3; RO = Yey 
F6D6 CF 03; 50 § 904 1 2¢ POLYF RE TUR as INDLN, SINDTB ; RO = PCY*Y) 
38 54 44 SDA 1 28 MULF R4, R 3 RO = YeP(yry) 
54 2 00 ef A 90D 1 SUBW #71300, R4 3 R4 = ¥22*-6 
0 4 eb ae 1 28 ADOF R4, R ; RO = SIN(Y) 
0 9E 1 RETURN: RSB 
9E6 1860 
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-SBTTL DEGENERATE SOLUTIONS 


P_ONE 
50 08 30 MOVF #1, RO ; Answer is 1 
RSB 
N_ONE: 
50 0000C080 8F 30 2 MOVF #-1, RO ; Answer is -1 
0 RSB 


NOUS WN "OO OONOAUE Wy 


; Underflow; if user has FU set, signal error. Always return 0.0 
MOVPSL 
CALLS 


WOOOOOOODOOOOOOOOOOOOO 
sel al etl esl esl eles esl esl esl asl eal 


DW SPMNNININININ— DB PS FP OPAAAAAD 


: = = user's or jacket routine's PSL 
= TRUE it JSB from jacket routine 

; branch if user ote ag 

et user PSL saved a7 fA 

0 = result. LIBSSIGNAL mi save in 


R 

0, G*MTHSSJACKET_TST 
BLEC =O sR 

MOVZWL SF$W_SAVE_PSW(FP), R2 
CLRL RO 


COooCoooo 


04 
52.04 


6 
7 
7 
7 
7 
7 
76 UNFL: 
7 
7 
7 
52 8 
00000000 ' GF 9? 
AD 
50 : 
i 
8 
9 


AO 10$: : 
aN H : mikirr MCH_ RO/R1 so any handler can 
0D 52 06 €1 QA04 : BBC #6, R2, 208 has user enabled fleeting ¥ under f low? 
DD OA08 PUSHL (SP) ; yes, return PC from special routine 
7E 00°8F 9A QAOA 4 MOVZBL #MTHSK_FLOUNDMAT, -(SP) | ieee code for hardware floa 


tin 
erflow convert tc Mite FLOUNDMAT 
: (32-bit VAX-11 oucapt ten Code) 
; signal (condition, P 
; return 


00000000 ° GF 


02 


1891 
B OA0 1833 CALLS #2, G*MTHSS$SIGNAL 
5 18935 20$: RSB 

1894 
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~ END 
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EVAL 


UVXSSINCOS 
Symbol table 


SINLENR 


SMA IN 
SMALL-SINCOS 
SMALL -SINCOSD 
SMALL-SIND 


UAB JUST 
W_MAX_WEIGHT 
W_TERA_WEIGHT 


X_1_0V_45 


PSECT name 


37124 bytes 


Initialization 
yy processing 


s 1 
Syabol table sort 


as table output 
Psect synopsis output 
Cross-reference output 
Assembler run totals 


Limit was 1500 pages 
pages) of virtual memory were used to buffer the intermediate code. 
There were 20 pages of symbol table space aceqcates to hold 196 non-local and 49 local symbols. 
1955 source Lines were read in Pass 1, §2 2. 
10 pages of virtual memory were used to define 9 macros. 


The working sst 


Macro library name 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 


Page faults 


. 7% 
; Floating Point Sine, Cosine and Sincos .« 


~138¢ OG: 96: 1 ae Se v04-00 Page 
MTHRTL.S JUvx XSINCOS.MAR;1 
a 
CR 
eR 
9R 
see 
see 
; 96 R 0 
Do009ES RO 
sae 
0904900 
00001000 
00000004 
00000086 
See esc oe eas eoeeeas + 
: Psect synopsis : 
Allocation PSECT No. Attributes 
00000000 ¢ 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
00000000 < Q.) O01 ¢ 1.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
QOO00A16 ( 2582.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


CPU Time Elapsed Time 
36 09:00:00.18 00:00:00.47 
126 00:00:00.69 00:00:04.65 
211 00:00:05.61 00:00: 14.68 
0 89:82:00.6 Bo? 0200 - 6 
327 00:00:03.96 00:00:1 06 
17 00:00:00.15 00:00:00.9 
2 00:00:00.02 00:00:00.02 
0 BF 808-88 88:80:00 -08 
721 00:00:10.80 00:00:33.19 


producing 34 object records in Pass 


ED 


3 Floating Point Sine, Cosine and Sthees Vo-see Pe -1388 99: 96: 12 PAr/vas Be ey ver =) B nan: Page 8, 


Hh tha ees Run Statistics 
131 GETS were required to define 5 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:UVXSINCOS/OBJ=OBJ$:UVXSINCOS MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) +MS 


NT CORPORATION 
PROPRIETARY 
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